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I.  PHYSICAL  APPROACH  TO  THE  PROBLEM 

There  exists  a  wealth  of  experimental  evidence  establishing 
the  proposition  that  irradiation  with  electromagnetic  (EM)  fields 
affects  gene  expression.  The  research  described  here  represents 
the  main  features  of  a  program  designed  to  develop  an  understand¬ 
ing  of  this  at  two  levels: 

(1)  Phenomenologically  we  have  developed  what  we  term  a  "deple¬ 
tion  model"  to  account  for  the  gross  features  of  the 
dynamical  effect  of  EM  fields;  and 

(2)  At  a  molecular  level  we  have  designed  and  are  performing  a 
set  of  experiments  aimed  at  elucidating  the  physical 
character  of  the  interaction  mechanisms  by  which  EM  energy 
is  coupled  into  DNA  molecules. 

In  the  paragraphs  that  follow  we  report  on  the  accomplishments 
and  progress  in  these  two  areas.  We  consider  first  the  dynamical 
modeling  efforts  and  then  describe  the  spectroscopic  work  than 
that  has  been  and  is  being  carried  out.  We  also  describe  the 
design  and  engineering  of  new  systems  for  achieving  desired 
irradiation  conditions  of  biological  systems. 

1 .  Dynamical  Modeling 

A.  The  Depletion  Model 

A  number  of  experimenters  (ourselves  included)  have  reported 
that  the  exposure  of  cells  to  relatively  low  intensity  pulsed  and 
low  frequency  electromagnetic  fields  results  in  a  transient 
augmentation  of  transcriptional  activity.  Surprisingly  the 

increases  have  been  seen  under  certain  irradiation  conditions  to 
exhibit  maxima  when  regarded  as  a  function  of  the  strength  of  the 
electromagnetic  fields.  We  have  developed  a  simple  multi-step 
chemical  reaction  model  that  accounts  for  the  principal  features 
that  are  observed  both  in  the  time  and  power  variation  of  the 
transcriptional  effects.  The  crucial  hypothesis  of  the  model  is 
the  supposition  that  the  direct  effect  of  cell  exposure  to 
electromagnetic  fields  is  an  increase  in  the  rate  constant 
characterizing  one  of  the  sequential  reactions  in  the  synthesis 
of  mRNA. 

We  assume  that  the  synthesis  of  messenger  RNA  and  the 
subsequent  protein  production  can  be  schematically  represented  by 
a  set  of  sequential  reactions  as 


A 


ki 


> 


X 


- > 

k2 


Y 
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Where,  for  example,  A  represents  the  nucleotide  pool  in  the  cell, 
and  the  process  A  -*  X  is  the  diffusion  controlled  migration  of 
these  nucleotides  to  the  neighborhood  of  the  DNA.  The  other 
steps  in  the  symbolic  reaction  scheme  represent  the  transcription 
process  (resulting  in  the  production  of  messenger  RNA,  denoted  by 
Y) ,  and  the  translation  process  (producing  proteins,  Z)  .  The  k's 
are  the  effective  rate  constants  for  the  various  steps. 

We  designate  the  concentrations  of  the  various  products  by 
x,  y  and  z  respectively  and  write  first  order  differential 
equations  governing  their  production: 

dx/dt  *  ki  A  -  k2X  (la) 

dy/dy  =  k2X  -  k3  y  (lb) 

etc . 

We  have  assumed  that  the  nucleotide  pool  A  constitute  a  reservoir 
that  is  rapidly  replenished,  and  also  that  the  back  reactions, 
e.g.,  Y  —  X,  can  be  ignored.  The  steady-state  situation,  is  des¬ 
cribed  by  dx/dt  =  dy/dt  =  •  •  •  =  0  and  x  =  xo ,  y  =  yo ,  ...  ,  from 
which  we  can  obtain 

xo  *  ki  A/k2  and  yo  =  k2  xo /k3  *  ki  A/ks  .  (2) 

Suppose  now  that  at  some  instant,  say  at  time  =  0  (when  steady- 
state  conditions  are  extant) ,  the  system  is  subjected  to  electro¬ 
magnetic  irradiation. 

We  assume  that  the  effect  of  the  electromagnetic  field 
is  to  alter  the  rate  constant  kx  ,  increasing  it  to  some 
new  value  kx  * . 

For  t>0,  the  concentrations  x,  y,  z,  ...  are  time  dependent.  It 
is  a  straightforward  matter  to  solve  the  differential  equations 
to  obtain  the  results: 


x  ( t ) 

and 


ktA 


[1  +  ( 


(3a) 


y  (t)  = 


k3-k2 


(k7}  [e"K; 


-  e'*3'] 


where  Ak  =  k2*  -  k2. 

These  two  equations  are  two  of  the  important  conclusions  of 
the  model.  Observe  from  Eq.  (3a)  that  x,  the  concentration  of 
nucleotides  near  the  DNA,  decreases  from  its  zero-field  value  of 
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ki A/ka  to  its  steady  state  value  ki A/k2 *  as  a  result  of  the 
increased  transcription  rate.  The  concentration  of  messenger 
RNA,  x(t) ,  first  rises  as  a  result  of  the  enhancement  of  the  rate 
ka ,  but  then  decays  eventually  to  its  basal  value  as  a  result  of 
the  depletion  of  the  nucleotides.  The  excess  mRNA  activity  above 
the  basal  value,  that  is,  y(t)  -  yo is  is  plotted  in  Figure  1  for 
several  values  of  the  parameter  A k/k2  (the  values  for  the  other 
parameters  are  given  in  Table  1) .  We  assume  that  Ak/k-  increases 
as  the  irradiating  field  strength  is  increased. 


TABLE  1.  Parameters 

used  in  constructing 

Figs.  1-5. 

Figs.  1-3 

Figs.  4-5 

A  (arbitrary 
units) 

10 

10 

ki  (min* 1  ) 

1.0 

1,0 

k2  (min* 1 ) 

0.001 

0.01 

k2  *  (min*  1  ) 

0.008 

0.040 

0.100 

0.1 

k3  (min*  1  ) 

0.800 

0.2  I 

Observe  in  Fig.  1  that  for  an  irradiation  time  of  20  minutes 
the  response  (i.e.,  y,  the  mRNA  concentration)  is  largest  for  the 
intermediate  exposure  conditions.  If  the  only  "data"  that  were 
presented  were  these  30-minute  exposure  results  considered  as  a 
function  of  irradiating  field  strength  (see  Fig.  2) ,  then  it 
would  be  quite  natural  to  infer  that  the  biological  system 
exhibits  an  enhanced  sensitivity  to  fields  near  the  intermediate 
value.  As  is  evident,  however,  such  a  conclusion  would  make 
little  sense  in  view  of  the  time  evolution  of  the  response  as 
shown  in  Fig.  1.  It  is  tempting  to  speculate  that  many  of 
the  so-called  "power  windows"  that  have  been  reported  in  the 
literature  might  be  explained  in  a  fashion  similar  to  this. 


Observe  that  even  under  the  assumption  that  4k  a  E,  there 
are  no  obvious  features  in  the  response  curves  of  Fig.  1  which 
are  linearly  related  to  the  field  strength  (the  height  of  the 
peak,  for  instance  does  not  quintuple  in  going  from  the  inter¬ 
mediate  to  the  largest  field  strength).  However,  by  differen¬ 
tiating  Eq.  (3b)  we  can  obtain  the  result 


k,A 

kz 


ak 
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which  indicates  that  the  initial  slope  of  the  response  curves  in 
Fig.  1  is  directly  proportional  to  Ak,  and  thus  (by  assumption) 
to  the  electric  field  strength.  This  behavior  is  evident  in  Fig. 
3  in  which  the  early  region  of  Fig.  1  is  shown  in  detail. 


We  may  also  ask  about  how  the  system  responds  when  the  field 
is  switched  off.  If  we  hypothesize  that  the  immediate  effect  of 
this  is  to  cause  the  return  of  the  rate  constant  to  its  original 
value,  i.e.,  k,*  -*  k.,  then  we  can  obtain  the  desired  results 
simply  by  resolving  the  differential  equation  subject  to  the 
appropriate  initial  conditions.  Assume  that  following  the 
initial  imposition  of  the  field,  sufficient  time  passed  to  allow 
both  x  and  y  to  reach  their  steady  state  values.  At  this 
instant,  say  time  =  t* ,  the  field  is  suddenly  switched  off.  The 
time  dependence  of  x  and  y  for  tst*  is  then  given  by 


x(t)  -  [1  -  (7^) 

z  K2 


and 


(4a) 


y  (t) 


ktA 

*3 


k,A 

k3-k2 


[e 


e  ] 


(4b) 


where  At  =  t  -  t*.  The  complete  time  dependence  of  the  mRNA 
concentration  is  plotted  in  Fig.  4.  Note  particularly  that  the 
mRNA  concentration  dips  below  its  steady-state  basal  value 
following  the  removal  of  the  radiation  field,  reaches  a  minimum 
value,  and  ultimately  returns  to  its  basal  value  as  t  -*  od.  This 
"rebound  "  effect  has  not  as  yet  been  reported  in  the  literature. 
We  have  initiated  efforts  (now  under  way)  to  observe  this 
behavior . 


One  last  point  should  be  discussed.  We  suppose  that  the 
proteins  produced  (Z  in  the  reaction  chain)  in  the  translation 
process  are  relatively  long-lived,  in  the  sense  that,  their  own 
denaturing  proceeds  much  more  slowly  than  all  other  relevant 
dynamical  steps.  Then  it  is  a  fairly  simple  matter  to  obtain  the 
protein  concentration  z(t)  from  dz/dt  =  k3 y  as 

z  (t)  =  k3  /  dt '  y  (t ' )  ; 

0 


z(t)  =  kjAt  + 


ktA 

k,-k. 


’  Ak\  r  1^3. ,  4 _ ~y'z  1 


(£*)  [  «  ( l~e  2  )  -  ( 1-e  3  )]  . 

Kz  k„ 


(5) 


The  linearly  increasing  first  term  which  describes  the  normal 
production  of  protein  in  the  absence  of  external  irradiation.  If 
we  focus  on  the  second  term  (call  it  6z)  which  describes  the 
effect  of  the  irradiation,  it  is  clear  that  this  term  initially 
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rises  quadratically ,  reaches  an  inflection  point,  and  ultimately 
levels  off  at  a  limiting  value  of 

6z(t— oo)  =  5z„  =  k.A  (j^-  -  r~r )  • 

2  •»  2 

Observe  that  for  kz  *  >>  kz  ,  this  becomes  5z-  *  ki A/kz ;  the 

limiting  value  is  independent  of  the  strength  of  the  irradiating 
field,  i.e.,  of  the  change  in  k* .  The  variation  of  5z  is  shown 
in  Pig.  5.  This  result  means  that  one  cannot  linearly  extra¬ 
polate  large  dose  observations  to  lower  doses  (as  is  often  done) 
to  establish  "safe"  exposure  levels.  The  dependence  of  the 
response  (5z)  on  the  field  strength  is  just  much  too  complicated 
for  a  naive  linear  extrapolation  to  be  meaningful. 
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EFFECT  OF  EM  FELDS  ON  RNA  SYNTHESIS 


TBC  DURATION  OF  EXPOSURE  MNLTTES 


1. 


Excess  oRNA  activity  plotted  as  a  function  of  tine  [y(t)  vs.  t] 
following  the  switching  on  of  an  incident  electromagnetic  field 
computed  from  Eq.  (3b).  The  parameters  used  are  given  in  Table 
1,  with  the  curve  wjth  the  highest  peak  corresponding  to  the 
largest  value  of  k,  . 


EFFECT  OF  EM  FELDS  ON  RNA  SYNTHESIS 


EXPOSURE  TICOO  fcMJTES 


EXCESS  PROTEIN  PROOUCED 


AS  FUCnOU  OP  QECTRC  TUD tTKNQTM 


ELECTRIC  HELD  STRENGTH 


rigure  5.  The  dependence  of  Che  excess  procein  production  6 2®  as  a  function 
of  field  strength.  The  curve  Is  computed  from  Equation  (6)  under 
the  assumption  that  Ak  *  E. 
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B.  Molecular  Dynamics 

The  molecular  dynamics  computer  simulation  studies  were 
designed  to  provide  some  insight  about  a  number  of  features  of 
polyelectrolyte  solutions.  These  include  the  following: 

(1)  We  have  considered  the  energy  transfer  processes  that  occur 
in  the  interaction  between  a  charged  polymeric  molecule  in 
an  electrolytic  solvent  and  a  suddenly  impressed  electric 
field.  In  particular,  we  have  attempted  to  examine  the  non¬ 
resonant  excitation  of  internal  vibrations  of  the  polymer, 
resulting  from  collisions  with  ions  set  into  motion  by 
externally  impressed  fields. 

(2)  We  have  attempted  to  understand  the  distribution  (cluster¬ 
ing)  of  counter-ions  and  co-ions  around  a  single  electri¬ 
cally  charged  polymer  molecule  in  solution,  both  in  the 
presence  and  in  the  absence  of  externally  applied  electric 
fields.  We  have  examined  the  time-dependent  response  of  the 
ionic  distribution  to  the  switching  on  of  the  field. 

The  model  system  under  examination  was  composed  of  1000  atoms, 
six  of  which  were  harmonically  bound  to  form  a  linear  molecule. 
These  six  atoms  were  each  given  a  negative  electric  charge.  The 
molecule  is  immersed  in  a  liquid  solvent  of  994  atoms  of  which  50 
are  positively  charged  and  a  number  are  negatively  charged,  that 
number  being  chosen  to  ensure  charge  neutrality  of  the  system  (it 
was  a  variable  as  various  values  were  examined  for  the  charge  on 
each  polymer  atom) .  All  atoms  in  the  system  interacted  with  one 
another  via  Lennard-Jones  6-12  forces,  as  well  as  the  appropriate 
electrical  forces  for  the  charged  particles.  At  various  times  in 
the  modeling,  one  or  more  of  these  electrical  forces  were  turned 
off  in  an  effort  to  understand  the  individual  effects  of  the 
various  interactions. 

In  the  early  studies,  no  restrictions  were  placed  on  the 
angle  between  adjacent  bonds  in  the  polymer  molecule.  The 
individual  atoms  thus  enjoyed  extraordinary  freedom  of  movement. 
In  the  presence  of  the  positively  charged  counter  ions,  the 
molecule  tended  to  wrap  itself  rather  tightly  around  one  of  the 
ions,  and  then  to  attract  a  roughly  spherical  sheath  of  addition¬ 
al  positively  charged  ions  around  it.  In  this  configuration  the 
polymer  was  quite  effectively  shielded  from  further  interactions 
with  moving  ions,  the  whole  assembly  simply  moving  electrophoret- 
ically  in  the  impressed  field. 


To  stiffen  the  polymer  backbone  so  that  the  polymer  will 
couple  more  effectively  to  ionic  motions  in  the  surrounding 
media,  we  added  a  set  of  constraining  forces  to  the  molecule  such 
that  the  angle  between  adjacent  bonds  is  maintained.  This 
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accelerated  the  exchange  of  energy  between  the  internal  vibra¬ 
tions  of  the  polymer  molecule  and  the  field-induced  dynamics  of 
the  surrounding  medium. 

For  this  system,  we  obtained  the  surprising  observation  that 
for  field  strengths  of  2  and  4  (in  reduced  LJ  units)  the  energy 
associated  with  internal  vibrations  decreases  when  the  field  is 
switched  on,  indicating  a  coupling  of  energy  from  these  vibra¬ 
tions  to  the  molecule's  ionic  environment.  That  is,  energy  is 
transferred  out  of  the  molecular  vibrations  into  the  energy 
reservoir  of  the  fluid  as  a  whole.  While  this  leads  to  the 
possibility  of  athermal  non-resonant  absorption  of  energy  by  the 
system  as  a  whole,  it  does  not  imply  excitation  of  the  elec¬ 
trically  charged  polymer  molecules. 

When  the  MD  experiment  was  repeated  from  a  different  initial 
configuration,  but  under  the  same  temperature,  pressure  and 
density  conditions,  essentially  the  same  generic  result  was 
obtained.  That  is  again  the  energy  was  observed  to  couple  out  of 
the  vibrational  degrees  of  freedom  of  the  molecule  and  into  the 
translational  energy  of  the  solvent.  Though  there  were  some 
quantitative  differences  in  the  runs,  the  qualitative  conclusions 
were  unchanged. 

In  addition  MD  experiments  in  solutions  in  which  no  external 
temperature  control  is  imposed  were  made.  In  contrast  with  the 
earlier  results  which  were  obtained  under  quasi-isothermal 
conditions,  i.e.,  at  each  time-step  in  the  computation,  the 
velocities  of  the  solvent  molecules  and  ions  were  scaled  to 
maintain  their  total  kinetic  energy  as  a  constant.  In  the  MD 
experiments  reported  here,  the  polymers  in  solution  were  made 
rather  rigid  with  only  internal  vibrational  modes  being  allowed. 
It  was  felt  that  this  would  accelerate  energy  transfers  between 
the  molecule  and  its  environment,  an  important  consideration, 
since  significant  temperature  rises  are  expected  to  occur,  and 
once  this  occurs,  interpretation  of  the  results  is  difficult. 

Immediately  after  switching  on  the  field  the  results  in 
terms  of  energy  coupling  were  essentially  as  before — a  slow 
transfer  of  energy  out  of  the  molecular  degrees  of  freedom  into 
the  fluid  as  a  whole  was  observed.  As  before  the  switching  on  of 
the  field  effected  a  tighter  packing  of  counter  ions  around  the 
polyions,  but,  since  these  molecules  are  in  this  study  quite 
rigid,  no  significant  changes  in  their  structure  occurred.  At 
somewhat  longer  times  after  the  field  was  turned  on  the  tempera¬ 
ture  of  the  system  rose  significantly  and  eventual  coupling  of 
energy  into  the  intra-molecular  vibrations  was  observed. 
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Spectroscopic  Studies 


A.  Coupling  of  Microwave  Energy  into  the  Environment  of  DNA 

Molecules 

Experimental  results1-17  have  shown  that  the  dielectric 
spectra  of  polyelectrolyte  solutions  display  several  dispersion 
regions  corresponding  to  different  mechanisms  of  polarization 
relaxation.  In  various  theoretical  approaches,  such  phenomena  as 
electrode  polarization,  fluctuations  in  the  distribution  of 
counterions,  restricted  motion  of  "bound”  water  molecules,  among 
others  have  been  considered  as  candidates  to  explain  the  observed 
relaxation  processes.  In  the  paragraphs  that  follow,  we  review 
some  of  these  mechanisms  to  provide  context  and  background  for 
the  results  and  their  interpretation  that  are  reported. 

Theory  of  counterion  condensation 

Most  of  the  models  that  have  been  advanced  are  independent 
of  developments  in  the  theory  of  polyelectrolytes,  the  fundamen¬ 
tals  of  which  were  proposed  by  Manning16-17.  Of  primary  impor¬ 
tance  in  this  treatment  is  the  phenomenon  of  counterion  condensa¬ 
tion,  which,  together  with  ionic  screening  interactions  at  long 
range  form  the  basis  for  this  general  theory. 

Counterion  condensation  refers  to  a  specific  mode  of  binding 
of  counterions  to  a  polyelectrolyte  in  which  a  constant  relation¬ 
ship  exists  between  the  axial  charge  density  of  the  polyelectro¬ 
lyte  and  the  valence  of  counterions  which  is  valid  over  a  broad 
concentration  range  of  the  counterion  species.  The  concept  of 
counterion  condensation  is  markedly  different  from  association 
reactions  between  ion  pairs  occurring  in  dilute  small  electrolyte 
solutions,  in  which  the  extent  of  ion  binding  and  the  tightness 
of  the  complex  are  inseparable.  In  polyelectrolyte  solutions  in 
an  environment  of  a  single  counterion  species,  the  extent  and 
tightness  of  binding  are  uncorrelated.  The  extent  to  which 
condensed  counterions  are  dehydrated  or  localized,  that  is,  the 
tightness  of  binding,  can  be  expected  to  vary  both  with  the 
nature  of  the  charged  groups  on  the  polyelectrolyte  and  the 
counterion  species.  Condensed  counterions  which  are  fully 
dehydrated  are  tightly  bound  by  several  neighboring  charged 
groups  and  are  said  to  be  localized.  Condensed  counterions 
retaining  all  their  water  of  hydration  corresponding  to  the 
hydration  state  in  pure  water  are  said  to  be  delocalized. 
Condensed  counterions  with  hydration  states  between  the  two 
extremes  are  said  to  be  loosely  bound. 

Counterion  condensation  on  a  polyion  occurs  when  the  axial 
charge  density  is  greater  than  a  critical  value.  When  this 
condition  is  met,  condensation  continues  until  a  saturation  point 
is  reached.  At  this  point,  the  saturated  charge  fraction,  that 
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is  the  ratio  of  counterions  to  available  charged  sites,  is  less 
than  one  and  is  invariant  to  further  increases  in  the  ion 
concentration  in  the  bulk  electrolyte.  To  quantify  this  observa¬ 
tion,  a  dimensionless  structural  parameter  Q  is  defined  as 
proportional  to  the  charge  density, 

q  J/cb 

Q  -  -  ,  (1) 

kT 

where  q  is  the  proton  charge,  e  is  the  dielectric  constant  of  the 
bulk  medium,  k  is  Boltzmann's  constant,  T  is  the  absolute 
temperature,  and  b  is  the  average  axial  charge  spacing.  The 
value  of  the  saturated  charge  fraction  (qB*t/q)  depends  only  on 
the  valence  of  the  counterions  and  the  axial  charge  density  of 
the  polyelectrolyte  and  is  equal  to  (NQ)*1  ,  where  N  is  the 
valence  of  the  counterion.  This  relation  was  obtained  by 
minimization  of  the  total  free  energy  of  the  system  against  the 
parameter  0m  which  represents  the  number  of  associated  counter¬ 
ions  per  unit  charge. 

The  free  energy  of  the  system  has  two  contributions  which 
vary  with  9m :  g*  i  ,  given  by  the  total  work  of  charging  the 
polyelectrolyte  ionic  groups  against  their  mutual  electrostatic 
repulsions,  and  gBi«,  given  by  the  mixing  of  free  cations,  bound 
cations,  and  solvent  molecules  (assuming  the  polyion  is  negative¬ 
ly  charged).  By  itself,  g* i  would  lead  to  complete  binding,  that 
is  N0n  =1;  g.ix  alone  would  lead  to  ©n  =0.  The  equilibrium 
state  is  determined  from  minimization  of  the  sum  g«  i  +  g*u .  For 
polyelectrolytes,  however,  ge i  increases  without  bound  as  the 
concentration  of  simple  electrolyte,  Cm ,  in  the  bulk  goes  to 
zero.  Minimization  of  the  sum  gB i  +  gnx  leads  to  a  relation 
which  in  the  limit  of  Cm  vanishing  can  only  be  satisfied  by 
values  of  On  given  by 

NOn  =  1  -  (NQ)-1 .  (2) 

Since  the  normalized  effective  charge  is  defined  by, 


qB  •  t 

-  =  1  -  N0M  ,  (3) 

q 

where  qn*t  is  the  effective  charge  on  each  charged  group, 
equation  (2)  implies  that  the  charge  fraction  goes  as  1/NQ.  It 
can  be  shown  from  this  result,  that  for  Q  >  1  the  polyion  is 

surrounded  by  a  local  concentration  which  is  positive  even  as  the 
concentration  of  the  counterion  species  in  the  bulk  tends  to 
zero . 
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Dielectric  models  based  on  counterion  distribution  hypotheses 


Only  one  author  {Oosawa1 8  -  1 • )  has  suggested  using  Manning's 
theory  to  obtain  predictions  about  the  dielectric  behavior  of 
polyelectrolytes.  In  his  treatment,  the  low  frequency  dispersion 
was  explained  in  terms  of  the  relaxation  of  the  polarization  of 
condensed  counterions  induced  by  their  diffusion  along  the 
polyion.  The  high  frequency  relaxation  was  attributed  to 
polarization  perpendicular  to  the  polyion.  Oosawa  concluded  that 
the  dispersion  parameters  are  independent  of  the  length  of  the 
polyion  but  increase  with  concentration  suggesting  contributions 
of  loosely  bound  counterions. 

Jn  the  model  introduced  by  Van  der  Touw  and  Mandel,20  the 
polyelectrolyte  was  represented  as  a  sequence  of  charged  rod-like 
subunits  in  an  arbitrary  but  fixed  configuration.  The  mobility 
of  the  associated  counterions  along  a  subunit  was  assumed  to  be 
different  from  that  of  counterions  moving  between  subunits.  This 
model  predicts  the  occurrence  of  two  dispersions.  The  high 
frequency  dispersion  was  attributed  to  induced  dipole  moments  due 
to  the  distribution  of  counterions  along  the  subunits,  whereas 
the  low  frequency  dispersion  was  attributed  to  the  induced  dipole 
moment  of  the  entire  molecule  and  its  associated  counterions. 
Since  the  high  frequency  relaxation  is  determined  by  the  response 
of  ions  moving  within  each  subunit,  its  amplitude  is  independent 
of  the  length  of  the  macromolecule  and  consequently  of  its 
molecular  weight.  It  is  interesting  to  note  that  in  a  model  in 
which  the  polyion  was  represented  by  a  rigid  rod,  the  high 
frequency  dispersion  was  absent. 

Grosse21  (as  part  of  the  work  at  CUA  described  below,  and  in 
an  attempt  to  explain  data  obtain  in  our  laboratories)  considered 
a  theoretical  model  specifically  developed  for  the  case  of  DNA 
solutions,  but  also  applicable  to  other  polynucleotides.  His 
work  concentrated  on  explaining  the  nature  of  the  high  frequency 
relaxation.  The  DNA  solutions  were  represented  by  suspensions  of 
charged  cylindrical  rods  with  associated  "condensed"  counterions 
from  the  bulk  electrolyte.  The  response  of  the  "condensed" 
counterions  to  a  high  frequency  field  was  taken  to  be  equivalent 
to  that  of  a  highly  conducting  layer  surrounding  the  cylindrical 
particle.  With  this  model,  Grosse  calculated  the  dielectric 
properties  of  the  equivalent  homogeneous  particle  using  the 
mixture  formulas  derived  by  Weiner22  and  Rayleigh23.  For  prolate 
spheroids,  which  can  approximately  represent  long  cylindrical 
particles,  Fricke's  extension  of  the  Maxwell-Wagner  model  to 
elliptical  particles2 «  predicts  the  occurrence  of  two  relaxa¬ 
tions,  one  at  low  frequencies  and  one  at  high  frequencies.  Using 
the  general  method  of  Reynolds  and  Hugh2 5  Grosse  calculated  the 
parameters  of  the  latter  relaxation  showing  that  they  are 
independent  of  any  specific  absorption  mechanism  of  the  particle. 
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but  depend  on  the  fraction  of  ions  which  are  condensed  on  its 
surface. 


Dielectric  models  based  on  bound  water 

Experimental  studies  extending  beyond  100  MHz  have  pointed 
to  the  existence  of  another  relaxation  between  the  g  dispersion, 
and  the  relaxation  of  pure  water  (Schwan,11  Grant  et  al . 1 2 ~ 1 8 ) • 
This  relaxation,  called  the  5  relaxation,  which  has  been  observed 
in  proteins  such  as  hemoglobin  and  myoglobin  and  in  DNA  has  been 
attributed,  by  Grant  and  co-workers1 s  for  instance,  to  the  effect 
of  bound  water.  For  the  case  of  myoglobin,  Grant  pointed  out 
that  the  experimental  data  can  be  fit  to  obtain  two  Cole-Cole 
dispersions,  two  Debye  type  dispersions  and  one  Cole-Cole,  or 
four  Debye  type  dispersions.  With  additional  information  from  X- 
ray  diffraction  regarding  the  shape  and  size  of  the  molecule,  and 
using  mixture  relationships  based  on  the  Maxwell-Wagner  equations 
(discussed  later)  with  the  appropriate  model  for  the  hydrated 
molecule,  he  concluded  that  the  correct  representation  was  the 
one  yielding  four  Debye  type  dispersions.  The  first  two  disper¬ 
sions  correspond  in  this  model  to  the  a  and  g  relaxations,  the 
last  dispersion  was  forced  to  fit  that  of  pure  water,  and  the 
remaining  dispersion  which  occurs  in  the  low  GHz  range  was 
attributed  to  bound  water.  Grant  also  pointed  out  that  the  bound 
water  model  can  indeed  account  for  most  if  not  all  of  the 
observed  relaxations. 

The  Maxwell-Wagner  model  for  heterogeneous  mixtures  of  spherical 
particles 

If  a  polyelectrolyte  solution  is  viewed  as  a  heterogeneous 
mixture,  the  3  dispersion  can  be  attributed  to  interfacial 
polarization  also  referred  to  as  the  Maxwell-Wagner  effect- 26,27 
It  is  well  known  that  heterogeneous  systems  of  spherical  par¬ 
ticles  dispersed  in  a  continuous  phase  have  relaxations  due  to 
interfacial  polarization  occurring  at  the  boundaries  between  the 
components  of  the  mixture.  The  Maxwell-Wagner  derivation 
considers  a  system  containing  spheres  of  radius  a,  permittivity 
EP  and  conductivity  oP  sparsely  distributed  in  a  medium  of 
permittivity  £■  and  conductivity  oB .  The  system  is  acted  upon  by 
an  electric  field  E  *  Eo  exp(iwt) .  The  derivation  is  carried  out 
in  two  steps;  first  Laplace's  equation  is  solved  to  determine  the 
potential  outside  an  individual  particle  and  the  results  are 
extended  to  solve  for  the  potential  surrounding  a  sphere  of 
radius  R>>a  which  is  assumed  to  contain  all  the  particles,  and  is 
surrounded  by  the  continuous  phase.  Secondly,  the  same  problem 
is  solved  for  a  sphere  of  radius  R,  with  the  properties  of  the 
heterogeneous  system,  immersed  in  a  medium  with  the  properties  of 
the  continuous  phase.  The  solutions  obtained  for  the  potential 
outside  the  large  sphere  are  equated  to  obtain  for  the  specific 
admittance  K  of  the  mixture  (Van  Beek2 8 )  the  relation 
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K  -  Ka 


2Ka  +  Kp  +  2v  (  Kp  —  Ka  ) 
2  Ka  +  Kp  —  V  (  Kp  ~  Ka  ) 


(4) 


where  v  =  N (a/R) 3  is  the  volume  fraction  of  the  dispersed  com¬ 
ponent,  K  =  o  +  iwto £  is  the  specific  admittance  or  generalized 
complex  conductivity,  and  to  is  the  permittivity  of  vacuum.  This 
equation  can  also  be  expressed  in  terms  of  complex  permittivities 
e«  »  e  +  o/iw£o  as  shown  by  Grosse29  as  well  as  in  terms  of  other 
physical  parameters. 

The  frequency  dependence  of  the  permittivity  £  and  the 
conductivity  o  as  obtained  5  from  equation  (4)  shows  a  single 
relaxation  behavior.  Represented  in  a  Cole-Cole  plot  in  terms  of 
real  and  imaginary  parts,  £"  vs.  e'  and  o"  vs.  o'  appear  as 
semicircles.  Using  a  single  relaxation  Debye  model  with  e'  and 
e"  are  given  by 

£  8  “  £« 

£'  =  £-+  -  ,  (5) 

1  +  W  2  T  2 

and 

( £»  -  £-  )  WT 

t "  =  -  (6) 

1  +  W  2  T  2 

with  limiting  static  (e»)  and  infinite  frequency  (e->  p«r»it- 

t  1  V  1  t  1  «  S  , 

2  *  a  +  Op  +  2v  ( Op  —Ob  )  (  2  Oa  +Op  )  (  Ep  —  Eb)~  (2Eb+Ep)  (  Op  —Ob  ) 

£S  =£■  -  +3vob  - ; 

20b  +  Op  -  v(Op-Ob)  [20b  +  Op  -  V(  Op -Ob  )  ]  2 

(7) 

and 

2Eb  Ep  +  2v  (  £p  -  £■) 

E-  =  Ep  -  .  (8) 

2  Eb  +  Ep  -  v(Ep  -  Ea  } 

The  relaxation  time  t  is  given  by 

2£a  +  Ep  +  V  (  Ea  -  Ep  ) 

T  =  £o  -  .  (9) 

20b  +  Op  +  V  {  Ob  -  Op  ) 

In  this  case,  the  dielectric  parameters  are  independent  of 
particle  size.  It  must  be  noted  that  the  limiting  values  of  the 
permittivities  are  applicable  only  with  respect  to  the  Maxwell- 
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Wagner  effect.  When  other  dispersions  are  present  these  values 
represent  intermediate  constants  between  relaxations. 

Dintzis  et  al.«  postulated  that  for  solutions  of  polyelec¬ 
trolytes  cP=£m  and  Op >>o«  and  indicated  that  the  dielectric 
increment  ( £’-£»)  was  enhanced  as  the  difference  between  Op  and 
Ob  increased. 

Linear  polyions  can  be  approximated  as  prolate  ellipsoids  of 
revolution  with  large  values  of  the  axial  ratio.  O'Konski30 
considered  a  model  for  such  ellipsoids  surrounded  by  an  electri¬ 
cal  double  layer  characterized  by  a  surface  conductivity  assumed 
to  be  large  along  the  particle  surface  and  small  along  other 
directions.  He  derived  expressions  for  the  dielectric  increment 
and  the  relaxation  time  in  terms  of  components  of  the  depolariza¬ 
tion  factor  Aj  and  the  conductivity,  including  the  contribution 
of  the  surface  conductivity,  along  the  three  principal  axes  of 
the  ellipsoid: 

1  £j  / Ea  -  1  +  (oj/om  -1)  Aj 

£  —  -  V  £  a  I  -  (10) 

3  j  [1  +  (Oj  /Oa  "  1)  Aj  ]  » 

Ej  +  Ea  (Aj  “  1  “  1 ) 

T  =  £  o  -  (11) 

Oj  +  Oa  ( Aj  "  1  ~  1  ) 

In  this  case,  the  dielectric  properties  are  still  indepen¬ 
dent  of  the  size  of  the  particle,  but  depend  on  the  axial  ratio. 

O'Konski' s  model  was  not  adequate  to  explain  the  behavior  of 
systems  displaying  very  high  permittivities  at  low  frequencies. 
To  describe  these  cases,  Schwarz31  introduced  a  frequency 
dependent  surface  dielectric  constant  with  a  relaxation  deter¬ 
mined  by  the  diffusion  of  counterions  in  the  double  layer. 
Takashima3 2  extended  this  model  to  the  case  of  ellipsoids.  In 
both  cases,  the  Maxwell-Wagner  equations  still  hold  with  appro¬ 
priate  modifications.  Other  studies  have  followed  in  this  area; 
however,  since  they  concentrate  on  the  origin  of  the  low  fre¬ 
quency  dispersion  they  will  not  be  considered  here. 

( 1 )  Experimental  Research 

( a )  Methodology 

Various  types  of  instrumentation  have  been  used  to  measure 
dielectric  properties  of  solutions.  At  frequencies  below  100  MHz 
admittance  bridges  are  commonly  used  which  in  principle  allow  for 
measurements  down  to  dc.  In  this  case,  however,  there  are 
practical  limitations  particularly  with  conducting  solutions,  due 
to  electrode  polarization  effects,  which  generally  occur  at  low 
frequencies.  With  this  technique,  complex  permittivities  can  be 


20 


measured  by  placing  the  sample  between  two  plates  forming  a 
capacitor.  For  lossy  solutions  the  system  can  be  represented  by 
a  parallel  capacitance  CP  and  resistance  Rp .  The  complex 
admittance  of  the  circuit  is  given  by  Y*  ■  1/RP  +  iwCP .  The 
complex  permittivity  of  the  sample  e* ,  is  given  by 

t*  =  e'  -  i  (  £”  +  o/w£o  )  (12) 

where  e'  and  e"  are  respectively  the  real  and  the  imaginary  parts 
of  the  permittivity,  o  is  the  static  conductivity,  and  e0  is  the 
permittivity  of  vacuum.  The  complex  admittance  of  the  capacitor 
is  related  to  these  parameters  by  (Grant14) 

Y*  *  iwF  [ £  *  -  i(£"  +  o/weo)],  (13) 

where  F  is  a  geometrical  factor  for  the  capacitor.  The  param¬ 
eters  of  interest  £’,  and  em  can  be  calculated  from  CP  =  Fe‘  and 
1/RP  =  F(we"  +  a/to  ) .  As  can  be  seen  £"  can  not  be  directly 
obtained  except  for  cases  where  the  static  conductivity  is  zero. 
For  other  situations,  the  static  conductivity  must  be  indepen¬ 
dently  measured  for  instance  with  a  conductivity  meter. 

A  major  problem  for  measurements  at  lower  frequencies  is 
electrode  polarization.  This  effect  can  sometimes  be  a  problem 
even  in  the  low  MHz  range.  For  situations  where  electrode 
polarization  is  a  concern,  movable  electrodes  can  be  used  as 
suggested  by  Fricke  and  Curtis33  and  Schwan  and  Maczuck34- 

In  our  research  measurements  are  made  using  a  Hewlett 
Packard  4194A  Impedance/Gain-Phase  Analyzer  for  the  frequency 
range  between  100  Hz  and  40  MHz.  A  Hewlett  Packard  16451A 
dielectric  cell  with  movable  electrodes  is  used  in  this  range  to 
alleviate  the  electrode  polarization  problem.  Other  cell  con¬ 
figurations  are  occasionally  employed  when  necessary. 

Beyond  a  few  hundred  MHz  dielectric  measurements  can  no 
longer  be  made  with  the  impedance  bridge  technique.  Our  approach 
is  to  use  coaxial  lines  which  can  generally  be  utilized  between 
50  MHz  and  18  GHz.  This  technique  makes  use  of  the  reflected 
signal  from  a  sample  located  at  a  defined  reference  plane  along 
the  coaxial  line.  A  wide  range  of  sample  configurations  are 
possible.  In  all  cases  the  sample  is  placed  at  the  end  of  a 
coaxial  line  having  an  inner  diameter  a,  an  outer  diameter  b,  and 
a  characteristic  impedance  Zo .  The  sample  may  constitute  part  of 
the  coaxial  line  or  be  placed  at  the  end  of  the  line.  Two 
configurations  are  used  in  the  measurements:  a  modified  GR-900 
connector  (Stuchly39)  and  an  open  ended  coaxial  line.  The  GR-900 
cell  is  employed  up  to  2  GHz,  while  the  open  ended  coaxial  line 
is  employed  between  2  and  18  GHz.  A  Hewlett  Packard  8510  Network 
Analyzer  is  employed  to  determine  the  reflection  coefficient  in 
terms  of  the  scattering  parameters. 


21 


(b)  Areas  of  Study 


The  argument  regarding  the  interpretation  of  the  experimen¬ 
tal  results,  particularly  the  origin  of  the  high  frequency 
dispersion  still  goes  on  today.  One  objective  of  this  research 
is  to  use  the  fundamental  concepts  of  the  theory  of  polyelectro¬ 
lytes  proposed  by  Manning16-17  to  investigate  the  origin  of  the 
high  frequency  dispersion  of  specific  systems  of  synthetic  and 
biological  polyelectrolytes. 

Although  one  author  (Oosawa19)  previously  suggested  using 
Manning's  theory  to  describe  dielectric  behavior,  no  systematic 
investigation  has  been  carried  out  based  on  this  theory.  The 
parametric  nature  of  the  counterion  condensation  model  offers 
many  possibilities  to  design  such  a  study  to  test  this  model. 
The  fundamental  notion  is  that  the  localized  concentration  of 
ions  in  the  vicinity  of  the  polyion  chain  is  much  greater  than 
the  concentration  of  ions  in  the  bulk  elecLrolyte,  even  as  the 
latter  goes  to  zero.  Thus,  the  polyion  can  be  viewed  as  having 
two  types  of  counterions,  a  "condensed"  counterion  layer  whose 
concentration  is  (approximately)  invariant  with  the  concentration 
of  the  vbulk  electrolyte,  and  a  diffuse  counterion  layer.  The 
characteristics  of  these  regions  will  dictate  the  dielectric 
behavior  of  the  system.  There  are  several  parameters  which  may 
be  expected  to  influence  dielectric  behavior  as  predicted  by  this 
theory,  these  are:  the  bulk  electrolyte  species  and  therefore  its 
valence,  the  axial  charge  density  of  the  polyelectrolyte,  the 
nature  of  the  charged  groups  on  the  polyion,  and  consequently, 
the  level  of  hydration  of  the  counterions  and  the  tightness  of 
binding,  and  the  concentrations  of  polyelectrolyte  and  bulk 
electrolyte . 

Several  experimental  conditions  have  been  considered.  The 
effect  of  changes  in  the  level  of  hydration  of  condensed  counter¬ 
ions  are  investigated  by  testing  different  polyelectrolyte 
solutions  which,  while  having  the  same  concentration  of  condensed 
counterions,  exhibit  different  hydration  characteristics.  This 
study  (now  under  way)  will  yield  an  understanding  of  the  effect 
on  the  dielectric  behavior  of  the  tightness  of  binding  of 
counterions.  The  results  will  also  be  analyzed  in  terms  of  the 
bound  water  model  which  attributes  the  existence  of  the  5 
dispersion  to  hydration  layers  around  the  entire  macromolecule. 

Changes  in  the  permittivity  of  the  solvent  or  the  tempera¬ 
ture  affect  the  axial  charge  density  and  the  mobility  of  counter¬ 
ions.  The  effect  of  the  condensed  counterion  layer  on  the 
dielectric  behavior  is  under  investigation  by  increasing  the 
dielectric  constant  of  the  solvent  or  the  temperature  to  the 
point  where  no  counterion  condensation  occurs.  This  condition 
exists  when  the  thermal  energy  is  greater  than  the  electrostatic 
energy  of  interaction  between  charged  groups  on  the  polyion,  that 
is,  when  Q  S  1. 
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The  effect  on  the  nature  of  the  condensed  counterion  layer, 
and  consequently  the  dielectric  behavior,  of  deviations  from  the 
linear  conformation  of  the  polyion  is  being  studied  by  consider¬ 
ing  polymers  which  can  take  on  various  conformations  from  linear 
to  supercoiled.  For  the  latter  case,  the  Maxwell-Wagner  model 
for  spheres  is  expected  to  hold. 

It  has  been  suggested  by  Oosawa1 9  that  the  high  frequency 
dispersion  is  due  to  polarization  in  the  direction  perpendicular 
to  the  axis  of  the  polymer.  If  this  is  the  case,  the  dielectric 
increment  in  this  range  will  vary  with  the  diameter  of  the 
condensed  counterion  layer.  Experiments  to  test  this  hypothesis 
are  being  performed  by  varying  the  parameter  Q  which  specifies 
the  volume  occupied  by  the  condensed  counterion  layer,  and  thus 
the  effective  diameter  of  the  polymer. 

The  experimental  effort  is  concentrated  on  DNA  and  a  few 
other  representative  polymers  including  various  polynucleotides 
as  well  as  some  non-biological  polymers  such  as  polyvinyl 
sulfonate,  polyvinyl  sulfate,  and  polyphosphate. 
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(2)  Experimental  Results 


This  report  summarizes  the  results  of  the  second  year  of  a 
study  of  microwave  interactions  with  DNA  solutions.  The  initial 
interest  in  this  subject  was  prompted  by  reports  of  resonant 
absorption  of  microwaves  by  DNA1 * 2  .  It  was  reasoned  that 
relatively  small  amounts  of  energy  could  cause  significant  damage 
to  the  DNA  or  affect  its  behavior  if  the  resonances  were  undamped 
as  indicated  by  their  observation.  Recent  careful  measure¬ 
ments3  • 4  including  our  own3  have  shown  that  undamped  resonances 
are  not  present.  The  microwave  absorption  by  DNA  remains  of 
interest,  however,  because  of  the  fact  that  the  absorption  of 
energy  by  the  DNA  and  its  associated  counterions  may  be  propor¬ 
tionately  greater  than  in  the  surrounding  bulk  liquid.  The 
relatively  high  localized  absorption  may  affect  the  behavior  of 
the  DNA  and  its  ability  to  perform  its  functions. 

To  understand  the  microwave/DNA  interaction  we  have  found  it 
essential  to  characterize  thoroughly  the  chemical  composition  of 
our  DNA  solutions  with  particular  attention  to  ionic  conductors. 
We  have  extended  the  frequency  range  of  our  measurements  down  to 
the  kilohertz  range  to  be  able  to  completely  cover  the  relaxation 
processes  extending  into  the  microwave  region.  We  have  completed 
preliminary  measurements  of  the  dielectric  behavior  of  DNA 
solutions  from  a  few  kilohertz  to  18  gigahertz.  We  have  con¬ 
cluded  that  to  thoroughly  understand  the  observed  behavior  it 
will  be  necessary  to  perform  a  series  of  measurements  in  which 
the  ionic  species  and  concentration  in  the  DNA  buffer  is  varied. 
It  is  hoped  that  this  will  give  sufficient  information  to  allow 
us  to  choose  among  available  models  for  explaining  the  DNA 
dielectric  behavior.  This  knowledge  may  in  turn  lead  to  sugges¬ 
tions  of  how  electromagnetic  energy  may  affect  the  biochemical 
behavior  of  DNA. 

(a)  Microwave  Absorption  in  DNA  -  Measurement  Methods 

(i)  Automatic  Network  Analyzer  with  Open  Ended  Coaxial  Probe 

The  measurement  of  the  dielectric  properties  of  DNA  solu¬ 
tions  presents  several  formidable  obstacles  owing  to  several 
properties  of  DNA  which  include  its  relatively  low  solubility  and 
the  difficulties  involved  in  preparing  large  quantities  of 
purified  material.  The  first  of  these  difficulties  requires  that 
the  measurement  technique  be  highly  sensitive  and  the  second 
requires  that  measurements  be  made  in  relatively  small  volumes  on 
the  order  of  a  milliliter  or  less.  In  part  due  to  these  dif¬ 
ficulties  Edwards  et  al . 1  *  2  were  led  to  believe  that  they  had 
observed  resonant  absorption  in  DNA. 

To  measure  the  properties  of  DNA,  Edwards  e_t  al_.  used  and  we 
have  adopted  the  technique  of  reflection  measurements  using  an 
open  ended  coaxial  probe  and  an  automatic  network  analyzer  (ANA) . 
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The  details  of  the  measurement  technique  were  given  in  our  first 
annual  report8 .  In  order  to  observe  any  spectra  due  to  DNA  in 
dilute  solutions  (on  the  order  of  0.1%  by  weight)  it  was  neces¬ 
sary  to  subtract  a  large  background  due  to  the  buffer  solution 
and  observe  the  difference  spectrum  which  showed  a  resonance  like 
appearance.  In  the  frequency  range  where  the  resonances  were 
observed  any  geometrical  differences  between  the  reference  buffer 
measurement  and  the  sample  measurement  results  in  apparent 
resonances.  Such  geometrical  differences  could  be  a  small 
difference  in  sample  volume  or  a  minor  displacement  of  the 
measurement  probe  and  sample  container.  This  sensitivity  to  the 
geometry  of  the  measurement  is  an  unfortunate  consequence  of  the 
fact  that  to  extract  the  signal  of  interest  from  the  signal 
detected  by  the  network  analyzer  substantial  corrections  to  the 
signal  are  required.  These  corrections  are  required  to  elimi¬ 
nate  effects  of  signal  loss,  directivity  errors,  impedance 
mismatches  and  include  spurious  reflection  signals  such  as  from 
sample  containers,  connectors  etc.  Tf  these  errors  are  repro¬ 
ducible  then  they  may  be  corrected  but  if  they  change  from  sample 
to  sample  as  would  occur  if  sample  volume  is  changed  or  if  the 
probe  is  displaced  with  respect  to  the  sample  container  the 
corrections  will  be  in  error  and  a  noisy  resonant  like  appearance 
will  result  and  be  exaggerated  in  any  difference  spectrum. 

During  the  first  year  of  this  project  we  were  able  to 
duplicate  the  measurements  of  Edwards  et.  al.1>2  but  showed  that 
the  resonances  we  observed  were  simply  artifacts  which  could  be 
manipulated  by  geometrical  means.  Although  we  were  unable  to 
find  any  resonance  absorption  in  DNA  we  did  observe  what  appeared 
to  be  an  excess  absorption  which  was  largely  frequency  indepen¬ 
dent  below  at  least  9  GHz.  We  began  the  second  year  of  the 
project  with  the  intention  of  examining  and  explaining  this 
apparent  absorption  excess.  The  magnitude  of  the  excess  was  such 
that  if  it  were  associated  with  the  DNA  it  meant  that  a  dis¬ 
proportionate  fraction  of  electromagnetic  energy  falling  on  the 
samples  would  be  absorbed  by  the  DNA  and  its  associated  counter¬ 
ion  cloud. 

In  order  to  accomplish  this  task  it  was  necessary  to  improve 
on  the  measurement  technique.  We  had  to  eliminate  the  geometri¬ 
cal  effects  leading  to  the  apparent  resonances.  We  also  found 
that  it  was  necessary  to  control  the  temperature  of  the  probe  to 
prevent  gradual  drifts  in  the  signal  during  lengthy  measurements. 
Further  we  found  that  the  ionic  content  of  samples  did  not  agree 
with  the  expectations  of  the  suppliers  and  had  to  be  carefully 
determined  by  analysis  or  controlled  by  additional  sample 
purification  or  preparation  steps. 

To  demonstrate  that  spurious  reflections  from  the  sample 
container  could  be  eliminated  we  performed  calibration  and 
measurement  steps  in  increasingly  larger  containers  until 
movement  of  the  probe  in  the  sample  container  between  measure- 
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ments  no  longer  produced  resonances.  We  began  the  tests  in  our 
standard  1.5  ml  sample  tubes  and  found  that  to  completely 
eliminate  the  resonances  a  20  liter  container  was  required.  This 
would  not  be  practical  for  biological  samples. 

It  had  been  suggested  that  spurious  reflections  could  be 
eliminated  by  using  a  thin  walled  sample  container  surrounded  by 
an  absorbing  medium* .  We  experimented  with  a  number  of  con¬ 
tainers  including  finger  cots  and  condoms.  The  best  results  were 
obtained  with  the  thinnest  condoms  surrounded  by  a  saline 
solution  for  absorption.  With  such  an  arrangement  temperature 
control  and  handling  of  the  samples  presented  difficulties  and 
results  were  not  consistent  due  in  part  to  the  failure  of  such 
thin  walled  containers  to  maintain  any  fixed  shape.  We  concluded 
that  this  would  not  be  a  practical  approach. 

After  further  experimentation  we  found  that  if  great  care 
was  taken  it  was  possible  to  keep  the  measurement  geometry 
relatively  well  fixed  during  the  experiments  using  our  original 
sample  arrangement.  This  consisted  of  a  1.5  ml  polyethylene 
micro  centrifuge  tube  as  the  sample  container  which  is  placed  in 
a  concentric  temperature  controlled  cylindrical  jacket.  The 
sample  solutions  can  be  placed  in  the  sample  container  using  a 
syringe  fixed  with  a  small  plastic  tube.  The  syringe  can  also  be 
used  for  removing  the  sample  and  rinsing  the  container  without 
disturbing  either  the  container  or  the  sample  probe.  Using  such 
a  procedure  the  resonances  could  be  minimized  sufficiently  to 
make  useful  measurements.  In  Figs.  1  and  2  we  show  examples  of 
absorption  difference  spectra  made  using  the  old  procedures 
showing  resonances  and  the  new  procedures  showing  that  the 
resonances  have  been  largely  eliminated. 

A  second  source  of  error  in  our  early  measurements  using  the 
ANA  were  due  to  temperature  drifts  and  gradients  in  the  coaxial 
probe.  Even  after  lengthy  warmup  of  the  instrument  there  may  be 
a  temperature  gradient  between  the  end  of  the  coaxial  probe 
attached  to  the  instrument  and  the  end  immersed  in  the  sample. 
This  is  particularly  true  for  measurements  far  from  room  tempera¬ 
ture.  We  found  that  this  problem  could  be  minimized  by  thermally 
insulating  the  probe  from  the  point  of  connection  to  the  instru¬ 
ment  to  within  a  few  inches  of  the  sample  end.  Near  the  sample 
end  we  added  a  temperature  controlled  jacket  circulating  liquid 
at  the  same  temperature  as  the  sample  and  calibration  liquids. 
This  produced  a  stable,  reproducible  temperature  gradient  in  the 
probe  after  equilibrium  was  reached.  Equilibration  times  on  the 
order  of  60  minutes  were  typical. 


(ii)  Low  Frequency  Measurements-Impedance  Analyzer 

After  concluding  that  resonances  were  not  present  we  began 
to  concentrate  our  efforts  on  the  apparent  excess  conductivity  of 
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our  DNA  samples  and  on  the  known  but  not  fully  understood 
relaxation  spectra  of  DNA  solutions.  Relaxation  spectra  are 
generally  quite  broad,  extending  over  many  decades  in  frequency 
and  thus  it  was  necessary  to  extend  our  measurement  frequency 
below  the  45  MHz  available  with  the  network  analyzer.  We 
obtained  a  Hewlett  Packard  model  4194A  impedance  analyzer  (an 
automatic  impedance  bridge)  with  a  frequency  range  of  100  Hz  to 
40  MHz  allowing  us  to  cover  the  dielectric  spectrum  from  100  Hz 
to  26.5  GHz  with  two  instruments. 

To  make  conductivity  measurements  at  very  low  frequencies  we 
purchased  a  small  conductivity  probe  capable  of  making  measure¬ 
ments  with  sample  volumes  on  the  order  of  100  y liter.  This 
allowed  us  to  determine  a  conductivity  very  near  the  DC  conduc¬ 
tivity  of  the  sample  solutions.  The  probe  was  carefully  calibra¬ 
ted  using  saline  solutions  mixed  in  our  laboratories  and  for 
which  the  concentrations  were  verified  by  chemical  analysis.  We 
found  that  the  nominal  cell  constant  supplied  with  the  probe  by 
the  manufacturer  was  in  error  by  more  than  50%.  The  conduc¬ 
tivities  of  the  saline  solutions  were  calculated  using  an 
empirical  equation  due  to  Stogryn7 . 

To  make  dielectric  measurements  in  the  low  frequency  region 
we  constructed  a  cell  similar  in  design  to  one  used  by  Foster8  . 
This  cell  is  essentially  a  parallel  plate  capacitor  with  the 
sample  held  between  the  plates.  It  has  platinum  black  electrodes 
with  a  surface  area  of  0.7  cm2  and  for  which  the  spacing  can  be 
varied  from  0  to  2.0  cm.  The  capability  of  variable  electrode 
spacing  allows  one  to  make  various  measurement  corrections  for 
errors  such  as  those  due  to  electrode  polarization.  Further,  one 
can  within  limits,  by  changing  the  plate  separation,  adjust  the 
capacitance  and  resistance  of  the  cell  to  give  values  for  which 
the  impedance  analyzer  has  good  sensitivity.  The  cell  can  be 
temperature  controlled  by  circulating  temperature  controlled 
liquid  through  the  body  of  the  cell. 


(b)  Measurements  on  DNA 
(i)  Dielectric  relaxation 

For  samples  of  purified  plasmid  prepared  at  the  NIH  as  well 
as  commercially  purified  calf  thymus  DNA  we  made  measurements  of 
complex  dielectric  properties  from  a  few  kilohertz  to  18.5  GHz. 
DNA  concentrations  ranged  from  0.05%  to  0.172%  by  weight.  The 
measurements  above  45  MHz  performed  with  the  ANA  reproduced  those 
made  and  reported  previously  in  last  years  annual  report  with  the 
exception  that  the  formerly  observed  resonances  were  removed 
using  our  improved  measurement  techniques.  After  the  effects  of 
temperature  drifts  were  eliminated  and  with  the  addition  of  our 
low  frequency  bridge  measurements  we  were  able  to  verify  the 
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presence  of  a  relaxation  between  1  MHz  and  1  GHz  which  had  been 
observed  by  Takashima  ejt  al.9 

The  evidence  for  this  relaxation  can  be  seen  in  the  Fuoss- 
Kirkwood1 0  plot  in  Fig.  3.  £■ "  is  taken  to  be  the  value  of  the 
dielectric  loss  for  water  at  its  relaxation  frequency  and  e"  is 
the  dielectric  loss  for  the  DNA  solution.  In  this  plot  the 
approximately  linear  behavior  at  high  frequencies  is  due  to  the 
relaxation  due  to  water  near  17  GHz.  The  change  in  slope  at  the 
low  frequency  end  of  the  plot  is  an  indication  of  a  second  small 
relaxation  which  is  attributed  to  the  DNA.  To  understand  the 
nature  of  this  relaxation  we  believe  it  will  be  necessary  to 
perform  a  series  of  measurements  on  samples  as  a  function  of 
temperature,  DNA  concentration,  ionic  species  and  ionic  con¬ 
centration.  These  measurements  have  not  been  made  as  yet  and  we 
can  only  say  that  our  observations  are  consistent  with  those 
published  by  Takashima. 

The  importance  of  studying  the  effect  of  the  ionic  content 
of  the  solutions  is  that  we’  believe  that  the  dielectric  behavior 
of  the  DNA  is  intimately  affected  by  the  ion  cloud  surrounding 
the  DNA.  By  changing  the  ionic  concentration  and  valence  of 
these  ions  we  hope  to  obtain  sufficient  data  to  be  able  to 
explain  the  combined  effects  of  DNA,  its  surrounding  ions  and  any 
electromagnetic  radiation. 

(ii)  Conductivity  excess 

Our  early  measurements  of  microwave  absorption  in  our  DNA 
samples  indicated  an  excess-  absorption  due  to  a  conductivity 
excess  which  was  independent  of  frequency  below  at  least  9  GHz. 
The  magnitude  of  the  excess  varied  from  sample  to  sample  but  was 
sufficiently  large  that  it  was  thought  that  its  presence  could 
contribute  a  possible  damage  mechanism  for  DNA  due  to  localized 
absorption.  The  supplier  of  the  DNA  was  under  the  impression 
that  he  was  always  giving  us  samples  in  a  standard  buffer  having 
5  mM  NaCl,  0.1  mM  EDTA  and  lOmM  Tris-HCl.  Since  the  ionic 
content  of  the  samples  was  all  supposed  to  be  the  same,  we  at 
first  suspected  that  the  DNA  itself  in  the  solutions  was  causing 
the  conductivity  excess.  Grosse  proposed  a  model  in  which  the 
counterions  might  see  the  DNA  molecules  as  a  high  conductivity 
surface  with  the  net  result  being  an  increased  overall  solution 
conductivity.  Preliminary  calculations  based  on  this  model  using 
preliminary  experimental  parameters  were  not  able  to  explain  the 
magnitude  of  the  observed  conductivity  excess. 

Before  such  a  model  could  be  completely  tested  we  had  to 
account  for  all  possible  sources  of  normal  conductivity.  We 
first  attempted  to  verify  the  ionic  content  of  the  solutions 
using  chemical  analysis.  The  samples  were  analyzed  for  Na  and  K 
using  a  DC  plasma  instrument.  We  discovered  that  the  ionic 
content  of  the  solutions  was  in  general  several  times  that 
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specified  by  the  supplier  but  was  still  insufficient  to  explain 
the  observed  conductivity. 

The  supplier  suggested  that  triethylamine-acetate  (TEA) 
which  was  used  in  his  purification  procedures  might  contribute  to 
the  conductivity.  TEA  is  volatile  and  should  have  been  removed 
in  one  of  the  concentration  steps  but  there  was  a  possibility 
that  a  small  amount  might  remain.  We  measured  the  conductivity 
of  TEA  and  determined  the  concentrations  necessary  to  cause  the 
excess  conductivity.  Attempts  to  verify  or  disprove  the  presence 
of  TEA  in  the  samples  using  IR  absorption  were  unsuccessful  due 
to  the  small  signal  for  TEA  and  large  background  due  to  water  at 
the  possible  concentrations.  The  concentrations  of  TEA  required 
to  explain  all  of  the  conductivity  excess  in  the  samples  was, 
however,  greater  than  used  in  the  processing  of  the  samples  so  it 
was  concluded  that  TEA  could  not  be  a  major  source  of  the  excess. 

Toward  the  end  of  the  current  year  the  use  of  an  ICP  mass 
spectrometer  became  available  for  chemical  analysis.  Survey  runs 
on  the  samples  indicated  that  several  additional  ions,  Cs,  Ba,  Cu 
and  Mg,  in  the  samples  were  contributing  to  the  conductivity 
excess.  When  the  conductivity  of  the  known  ion  species  plus  that 
due  to  EDTA-tris  were  added  together  the  conductivity  was  still 
only  partially  explained.  Depending  on  the  sample,  the  remaining 
conductivity  excess  ranged  from  a  few  percent  to  more  than  50%. 
Discussions  with  our  chemical  analysis  group  convinced  us  that  it 
would  not  be  possible  to  accurately  identify  all  possible 
conductive  species,  especially  any  organic  ones.  Further,  it  is 
clear  from  examining  conductivity  models,  that  with  so  many 
possible  minor  species  contributing  to  conductivity  it  is  not 
possible  to  calculate  the  conductivity  with  great  confidence.  We 
determined  that  another  approach  was  needed  to  determine  the 
effect  if  any  of  DNA  on  the  conductivity.  We  decided  that  the 
best  approach  would  be  to  reduce  the  ionic  concentration  in  the 
DNA  solutions  as  much  as  possible  by  dialysis  and  then  add  back 
known  concentrations  of  NaCl. 

We  experimented  with  the  effectiveness  of  dialysis  in 
removing  unwanted  ions.  By  dialyzing  for  periods  up  to  one  week 
against  pu^e  water  we  found  that  we  were  able  to  reduce  the 
conductivity  of  the  samples  to  less  than  0.005  S/m  compared  to 
about  0.5  S/m  for  undialyzed  samples.  Chemical  analysis  indi¬ 
cated  that  Na  levels  were  less  than  0.5  ppm  (2  X  10-9  molar). 
Other  ions  detected  and  their  concentrations  were  Mg  0.1  ppm,  Cs 
1  ppm,  and  Ba  .2  ppm.  The  estimated  conductivity  from  these 
known  ions,  assuming  Cl  as  the  common  anion,  was  5  X  10" *  S/m  or 
about  10  %  of  what  was  measured.  Either  the  model  we  are  using 
is  grossly  inadequate  for  calculating  the  conductivity  or  their 
are  undetected  conductive  species.  We  suspect  the  later  but  can 
not  prove  it.  The  equivalent  NaCl  concentration  needed  to 
explain  the  excess  conductivity  is  on  the  order  of  10  ppm.  It  is 
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not  unreasonable  to  assume  that  some  undetected  organic  or 
inorganic  species  are  present  at  such  low  levels. 

The  normality  of  the  DNA  solutions  typical  of  our  samples  is 
about  3  mN.  Since  dialysis  was  able  to  reduce  the  NaCl  con¬ 
centrations  in  the  DNA  solutions  to  about  1%  of  this  we  are 
assured  that  the  DNA  has  been  effectively  stripped  of  NaCl.  If 
we  add  back  NaCl  to  a  DNA  solution  by  dialysis  against  a  saline 
solution,  we  expect  the  Na  ions  to  preferentially  become  counter¬ 
ions  to  the  DNA  for  NaCl  normalities  less  than  that  of  the  DNA. 
If  DNA  somehow  increased  the  effective  conductivity  of  associated 
counterions  then  upon  adding  NaCl  to  the  solution  we  should  see  a 
greater  conductivity  increase  than  when  adding  that  concentration 
of  NaCl  to  pure  water. 

We  performed  a  set  of  experiments  beginning  with  a  DNA 
solution  that  had  been  dialyzed  against  distilled  water  for  seven 
days.  As  explained  above  we  found  a  residual  conductivity  of 
0.005  S/m  with  about  10%  being  due  to  known  concentrations  of 
remaining  salts.  We  then  dialyzed  the  solution  against  saline 
solutions  of  0.0015,  0.005  and  0.015  normality  respectively. 

After  each  dialysis  the  conductivities  of  the  DNA  and  saline 
solutions  were  measured.  It  was  found  that  in  each  case  the 
conductivity  of  the  DNA  solution  increased  to  a  value  equal  to 
that  of  the  saline  solution  plus  its  original  unexplained  excess. 
Ignoring  that  small  excess,  the  conductivities  agreed  with  the 
value  calculated  from  the  empirical  formula  of  Stogryn  for  NaCl 
in  pure  water.  This  would  indicate  that  DNA  does  not  strongly 
affect  the  low  frequency  conductivity  of  associated  counterions. 
Since  counterions  are  expected  to  be  in  close  proximity  to  the 
DNA  the  local  salt  concentration  and  thus  the  localized  microwave 
energy  absorption  may  be  relatively  large.  This  could  be 
important  in  determining  the  effect  of  microwaves  on  the  behavior 
of  DNA.  To  further  understand  the  relationship  between  DNA  and 
the  counterions  it  is  necessary  to  study  the  dielectric  relaxa¬ 
tion  of  DNA  solutions  which  we  are  continuing  to  do. 

(c)  Conclusions 


The  results  of  our  dielectric  measurements  on  DNA  from  100 
Hz  to  18.5  GHz  do  not  indicate  any  unusual  behavior  for  the  DNA. 
There  is  evidence,  however,  for  dielectric  relaxation  between  1 
MHz  and  1  GHz.  Further  study  of  this  relaxation  is  desireable  in 
order  to  understand  the  process  of  electromagnetic  energy 
absorption  by  DNA  and  its  associated  counterions.  Since  the 
localized  energy  absorption  in  the  vicinity  of  the  DNA  may  be 
significantly  larger  than  in  the  bulk  it  is  possible  that  the 
biochemical  behavior  of  the  DNA  would  be  affected  in  some  way. 
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MICROWAVE  ABSORPTION 


FREQUENCY  («<Z) 


Figure  1.  The  absorption  difference  spectrum  for  two  water  samples.  The 
apparent  resonances  are  due  to  a  volume  difference  between  the  two  samples. 
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FREQUENCY  (OKZ) 


Figure  2.  The  absorption  difference  spectrum  for  two  water  samples  for 
which  the  sample  volumes  were  made  as  nearly  equal  as  practical.  The  resonance 
are  eliminated. 


DIALYZED  CALF  THYMUS  DNA 


7J  •  U  U  9J  M  tO  10u4 
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Figure  3.  FuosS-Kirkwood  Plot  of  a  calf  thymus  DNA  sample.  The  linear 
portion  of  the  curve  indicates  the  dominant  water  relaxation  near  17  GHz. 

The  change  in  curvature  at  low  frequencies  is  attributed  to  a  DNA  relaxation. 
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Theoretical  Efforts 


(a)  Localized  Microwave  Heating  of  DNA  Molecules  in  Solution 

The  conductivity  of  a  solution  of  DNA  is  determined  by  the 
free  ions  in  the  bulk  electrolyte  and  by  the  counterions  sur¬ 
rounding  the  charged  molecules.  The  total  numbers  of  these  two 
types  of  ions  can  be  of  the  same  order  of  magnitude,  so  that, 
when  and  external  field  is  applied  to  the  system,  the  generation 
of  heat  due  to  their  movements  can  also  be  comparable.  This 
raises  the  question  of  whether  a  localized  heating  of  the  DNA 
molecules  could  occur. 

If  such  an  effect  did  take  place,  it  wold  not  be  the  same 
for  all  the  DNA  molecules  in  the  solution  because  of  their 
different  orientation  with  respect  to  the  applied  field.  From 
the  biological  standpoint,  what  matters  is  not  the  average  value 
of  this  heating,  but  rather  the  maximum  temperature  increase 
expected  in  the  most  favorable  situation.  A  value  of  a  few 
degrees  would  then  indicate  a  possible  damage  mechanism  to 
biological  systems  exposed  to  electromagnetic  radiation. 
Speculation  about  such  localized  heating  effects  have  been  around 
for  some  time.1 

Calculations  of  the  microwave  heating  of  macroscopic  objects 
immersed  in  a  cooling  water  bath  have  also  been  previously 
reported.2  The  general  conclusion  is  that  a  differential 
temperature  rise  on  the  order  of  1°C  is  impossible  for  molecular 
size  objects.  Nevertheless,  this  conclusion  does  not  exclude,  a 
priori .  the  possibility  of  a  sizeable  heating  effect  in  DNA 
molecules  because  of  the  following  considerations: 

(a)  The  above  study2  was  analyzed  under  the  assumption  that  the 
cooling  fluid  in  contact  with  the  heated  object  had  a 
constant  temperature.  The  boundary  condition  is  well-suited 
for  experiments  in  which  a  circulating  fluid  is  used,  but 
greatly  overestimates  the  heat  transfer  in  the  case  of  a 
stagnant  medium. 

(b)  For  a  stagnant  bath,  the  static  problem  of  a  heated  object 

can  be  solved  for  simple  geometries.  For  a  sphere,  the 
temperature  increase  is  proportional  to  the  square  of  its 
radius,  and  it  has  a  negligible  value  for  molecular  size 
objects  and  any  reasonable  amount  of  heat  deposition.  The 
situation  is  totally  different  for  a  cylinder:  in  this  case 
a  static  solution  does  not  exist — the  temperature  of  the 
object  increases  continuously  with  time.  While  a  DNA 

molecule  cannot  be  strictly  represented  by  a  cylinder  of 
infinite  length,  so  that  it  must  eventually  attain  a  finite 
temperature,  this  steady-state  temperature  difference  should 
be  substantially  higher  than  that  of  a  sphere. 
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In  order  to  calculate  the  temperature  increase,  two  problems  must 
be  addressed: 


(a)  What  is  the  rate  at  which  energy  is  deposited  in  the  fluid 
as  a  function  of  the  strength  and  frequency  of  the  incident 
electromagnetic  field?  and 

(b)  What  is  the  efficiency  of  thermal  transfer  to  the  DNA 
molecule — and  thus  the  temperature  difference — as  a  function 
of  the  rate  of  heat  deposition? 

We  consider  a  linear  DNA  molecule  of  length  L  and  radius  R 
surrounded  by  a  layer  of  counterions  and  immersed  in  the  electro¬ 
lytic  solution.  The  dielectric  properties  of  the  molecule  are 
given  by  a  frequency  independent  permittivity  eP .  The  electric 
response  of  the  counterions  is  characterized  by  a  surface 
conductivity  r ,  while  the  electrolyte  is  characterized  by  a 
permittivity  £■  and  a  conductivity  o« . 

When  an  electric  field  Esw  exp  (iojt)  impinges  on  the 
system,  energy  is  dissipated  by  both  the  ions  from  the  bulk 
electrolyte  and  by  the  counterions.  The  power  deposited  on  the 
conductive  layer  is  determined  by  the  surface  conductivity  and  by 
the  local  field  strength,  which  strongly  depends  on  the  angular 
frequency  w. 

for  lower  frequencies,  the  highly  conductive  region  occupied 
by  the  counterions  tends  to  expel  the  field,  so  that  the  energy 
deposition  should  be  low.  Maximum  deposition  is  achieved  at  high 
frequencies  for  which  the  field  distribution  is  solely  determined 
by  the  permittivities  and  not  the  conductivities  of  the  com¬ 
ponents.  Since  we  are  interested  in  the  largest  possible  effect, 
we  shall  confine  out  attention  to  this  high  frequency  situation. 
Two  cases  are  to  be  considered,  those  of  molecules  oriented 
parallel  and  perpendicular  to  the  field. 

(a)  Molecule  Oriented  Parallel  to  the  Field 

In  this  case  (taking  L>>R)  the  molecule  does  not  perturb  the 
field,  so  that  the  field  amplitude  acting  on  the  counterions  is 
E«- ;  therefore  the  power  generated  per  unit  volume  of  the 
molecule  is 


Molecule  Oriented  Perpendicular  to  the  Field 

Here  the  field  inside  the  molecule  is  uniform  with  the  value 


tt  _  2€. 

E>  €P  +  c.  E- 
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Therefore  the  square  of  the  field  acting  on  the  counterions  is 


E.  (r=R) 


•e  +  e, 

P  a 


sin  0] 2  +  [g-^- 


cos  0] 


where  r  and  0  are  cylindrical  coordinates.  Integrating  the  power 
deposition  over  the  surface  of  the  cylinder  leads  to 


P 


(4) 


Since  typically  £a>>cP,  the  power  deposition  on  molecules 
oriented  perpendicular  to  the  field  is  about  twice  that  on 
molecules  aligned  with  the  field. 

The  maximum  expected  value  for  the  power  per  unit  volume 
that  can  be  generated  on  a  molecule  oriented  perpendicular  to  the 
field  can  be  estimated  as  follows.  The  amplitude  of  the  electric 
field  applied  to  a  system  is  limited  in  practice  by  the  electric 
breakdown  strength  of  air  (roughly  106  V/m) .  The  field  inside 
the  biological  system  is  generally  reduced  by  nearly  two  orders 
of  magnitude  due  to  the  high  permittivity  of  water.  This 
reduction  could  be  partly  compensated  by  the  "resonance  effect", 
which  occurs  when  the  specimen  size  is  comparable  to  the  wave¬ 
length.3  Thus  the  maximum  field  inside  the  sample  should  not 
exceed 


Em  -  (max)  a  103  V/m.  (5) 

The  surface  conductivity  f  can  be  estimated  from  the  number 
of  elementary  charges  per  unit  length  along  a  DNA  molecule4  (0.59 
•  1010  charges/m)  and  the  value  R  *  10- 9  m.  This  gives  a  surface 
charge  density  of  about  one  elementary  charge  per  100  A2  ,  which 
is  about  the  maximum  value  attained  by  surface  charge  densities 
in  solution.  Using  a  diffusion  coefficient  D  s  2  •  10-9  m2 /s 
then  leads  to 


r  a  5  •  10- •  S.  (6) 

The  heat  generated  per  unit  volume  in  a  DNA  molecule  could 
therefore  attain  a  maximum  value  of 


P  =  2  •  101  2  W/m3  .  (7) 

The  calculation  of  the  temperature  of  a  uniformly  heated 
object  in  a  stagnant  cooling  bath  can  only  be  carried  out 
analytically  for  a  few  simple  geometries.  For  a  sphere  the 
temperature  rises  monotonically  with  time  and  tends  to  an 
asymptotic  limit  which  coincides  with  the  solution  of  the  static 
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problem.  The  result  for  the  temperature  difference  between  the 
surface  of  the  object  and  the  initial  temperature  of  the  entire 
system  is3 


AT  =  PR2/3k  (8) 

where  k  is  the  thermal  conductivity  of  the  cooling  bath  (k  =  0.6 
W/m-deg  for  water) .  This  result  shows  that  for  spheres  with 
radii  less  than  one  1pm,  a  temperature  rise  of  1°C  cannot  be 
attained. 

For  a  spheroid,  the  final  temperature  difference  has  a 
simple  form  only  in  the  case  when  the  object  is  a  perfect 
conductor  so  that  the  temperature  insidef?  it  is  uniform.  In  that 
case 


(9) 

(10) 


This  result  diverges  when  e  tends  toward  unity,  showing  that  the 
final  temperature  difference  can  be  arbitrarily  large  for  a 
sufficiently  elongated  spheroid.  However  in  this  case  the  time 
that  is  required  to  reach  the  final  state  also  diverges. 
Consequently  it  is  necessary  to  examine  the  time  dependence  of 
the  temperature  change. 


The  result  takes  a  relatively  simple  form  if  was  consider  an 
infinitely  long  perfectly  conducting  cylinder.  The  asymptotic 
form  valid  for  t  >>  t  is 


AT(t)  =  ^  In  ^  (11) 

where  K  =  1.7811...  =  ey  (Y  =  0.57721  ...  =  Euler's  constant) 

and 


t  »  pCRVic,  (12) 

in  which  p  is  the  density  of  the  external  medium  and  C  is  its 
specific  heat.  For  DNA  molecules  in  an  aqueous  solution  r  » 
7  •  10- 12  seconds. 

Equation  (11)  clearly  constitutes  the  upper  limit  for  the 
temperature  increase  in  a  molecule,  since,  for  any  object  of 
finite  length,  the  temperature  must  reach  a  limiting  value.  The 
maximum  expected  temperature  difference  for  a  DNA  molecule  of 
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infinite  length  is  plotted  in  Figure  1  as  a  function  of  tine. 
The  two  horizontal  lines  correspond  to  the  limiting  values  for  a 
sphere  of  the  same  radius  and  for  an  ellipsoid  with  L/R  -  1000. 
The  conclusion  is  that  despite  the  large  amounts  of  power  per 
unit  volume  that  can  be  deposited  on  a  ONA  molecule  in  solution, 
no  appreciable  temperature  difference  can  ever  be  achieved 
because  of  the  rapid  transport  of  energy  away  from  an  object  of 
so  small  a  radius. 
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(b)  Microwave  Absorption  by  Suspensions  of  DNA-type  Particles 

A  theoretical  study  of  microwave  absorption  by  a  suspension 
of  charged  rod-like  particles  in  an  electrolytic  solution  has 
shown  that  due  to  the  concentration  of  ions  around  the  particles, 
the  absorption  of  the  suspension  is  frequency  dependent,  and  will 
always  exceed  at  high  frequencies  the  absorption  that  would  occur 
if  the  ions  were  uniformly  distributed. 

In  addition,  it  has  been  found  that  the  dielectric  behavior 
of  the  simplest  possible  model  of  a  suspension  of  charged 
cylindrical  particles  is  characterized  by  a  relaxation  in  the 
microwave  region.  The  parameters  of  this  relaxation  are  indepen¬ 
dent  of  any  specific  absorption  mechanism,  but  depend  on  the 
fraction  of  the  ions  which  are  condensed  on  the  surface  of  the 
particles . 

This  relaxation  has  its  origin  in  the  frequency  dependence 
of  the  field  distribution  inside  the  suspension.  At  low  frequen¬ 
cies  the  distribution  is  governed  by  the  relative  conductivities 
of  the  medium  and  the  particles.  The  field  tends  to  be  excluded 
from  the  counterion  layer  because  of  its  high  conductivity,  and 
is  concentrated  in  the  bulk  electrolyte.  This  reduces  the 
average  current  density  in  the  suspension,  so  that  the  conduc¬ 
tivity  falls  below  the  value  it  would  have  in  a  uniform  ionic 
solution. 

For  frequencies  above  the  relaxation,  the  current  densities 
are  unable  to  develop  appreciable  charge  densities  on  the 
boundaries.  The  field  distribution  is  then  independent  of  the 
conductivities,  depending  rather  on  the  permittivities  of  the  two 
phases.  The  field  is  enhanced  inside  the  particles  because  of 
their  relatively  lower  permittivity,  and  thus  the  field  near  the 
particle's  surface  (where  the  condensed  ions  are  located)  is  also 
enhanced . 

The  net  effect  is  that  the  energy  per  unit  volume  deposited 
near  the  surface  of  a  DNA  molecule  can  be  up  to  100  times  larger 
than  elsewhere  in  the  ionic  medium. 


(c)  Dielectric  Properties  of  Particles  in  Electrolytic  Solutions 

The  three  paragraphs  below  are  abstracts  of  completed 
theoretical  projects  which  are  pertinent  to  understanding  the 
mechanism (s)  of  microwave  absorption  by  solutions  of  DNA  or  other 
polyelectrolytes . 

i)  Permittivity  of  a  Suspension  of  Charged  Cylindrical 
Particles  in  an  Electrolytic  Solution  (by  C.  Grosse,  to 
be  published  in  the  Journal  of  Physical  Chemistry) 
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A  simple  model  recently  developed  for  the  low  frequency  permit¬ 
tivity  of  suspensions  of  charged  spherical  particles  in  electro¬ 
lyte  solutions  [C.  Grosse  and  K.  R.  Foster,  J.  Phys .  Chem.  91, 
3073  (1987)]  is  extended  to  the  case  of  cylindrical  particles 

oriented  perpendicular  to  the  applied  field.  Expressions  for  the 
conductivity  and  permittivity  as  functions  of  frequency  are 
deduced.  The  indicate  that  while  the  relaxation  occurs  at 
essentially  the  same  frequency  as  in  the  spherical  particle  case, 
the  low  frequency  behavior  is  totally  different.  The  low  fre¬ 
quency  conductivity  of  a  suspension  of  cylindrical  particles  has 
the  same  value  as  the  conductivity  of  the  electrolyte.  The  low 
frequency  permittivity  diverges  logarithmically,  while  its  value 
at  any  frequency  scales  with  the  conductivity  of  the  electrolyte. 

ii)  On  the  calculation  of  the  Dielectric  Properties  of 

Suspensions  (by  C.  Grosse,  to  be  published  in  the 

Journal  of  Ferroelectrics) 

The  dielectric  properties  of  suspensions  are  usually 
calculated  from  the  quasi-static  solution  of  the  Laplace  equation 
for  the  potential  around  a  single  particle  in  the  continuous 
phase.  The  results  are  only  straightforward  for  the  static 
conductivity  and  for  the  high  frequency  permittivity,  since  these 
expressions  only  involve  the  potentials  which  are  in  phase  with 
the  applied  field.  The  static  permittivity  and  the  high  fre¬ 
quency  conductivity  are  much  more  difficult  to  evaluate  since 
they  depend  strongly  on  the  limiting  behavior  of  the  out-of-phase 
terms,  which  actually  vanish  in  these  two  limits.  Using  argu¬ 
ments  based  on  the  stored  energy  and  on  the  dissipated  power,  it 
is  shown  how  expressions  for  these  properties  can  be  calculated 
sing  only  the  in-phase  terms  of  the  potential.  The  results  are 
much  easier  to  derive  than  those  using  the  classical  procedure, 
and  they  contribute  to  a  better  understanding  of  the  meaning  of 
the  relaxation  amplitudes.  A  further  advantage  of  this  method  is 
that  it  permits  one  to  obtain  the  limiting  solutions  for  the 
dielectric  properties  in  some  systems  (e.g.,  ellipsoidal  par¬ 
ticles  surrounded  by  a  layer  of  uniform  thickness)  for  which  an 
analytical  solution  of  the  Laplace  equation  does  not  exist  for 
all  frequencies. 


iii)  On  the  Extension  of  Maxwell's  Mixture  Formula  to 
Ellipsoidal  Particles  (by  C.  Grosse,  to  be  published  in 
the  Journal  of  Ferroelectrics) 

Fricke's  generalization  of  Maxwell's  mixture  formula  to  the 
case  of  ellipsoidal  particles  is  reexamined.  It  is  shown  that 
this  generalization  is  based  on  two  approximations  which  are  not 
justified  in  the  case  of  randomly  oriented  particles  with  high 
eccentricity.  A  new  expression  which  avoids  this  objection  is 
derived.  It  tends  to  Fricke's  result  at  low  volume  fractions, 
but  reduces  to  the  Maxwell  formula  for  very  high  concentrations. 
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When  particularized  to  the  case  of  needle-shaped  particles,  it 
shows  significant  deviations  from  Fricke's  formula,  which  are  in 
qualitative  agreement  with  the  extant  experimental  data. 
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B.  Coupling  of  Microwave  Energy  into  DNA  Molecules 

It  seems  clear  from  the  preceding  paragraphs  that  mechanisms 
exist  and  seem  to  be  operative  in  which  the  dissolution  of  a 
polyion  by  an  ionic  solvent  leads  to  enhanced  absorption  of 
electromagnetic  (microwave)  energy  by  the  solution.  The  question 
which  we  now  address  is  whether  microwave  energy  can  be  directly 
coupled  into  the  internal  degrees  of  freedom  of  such  polyions,  in 
particular  DNA.  In  the  microwave  frequency  region,  i.e.,  from  1 
to  about  20  GHz,  theoretical  considerations  lead  one  to  expect 
the  existence  of  longitudinal  vibrational  modes  of  the  DNA 
backbone.1’2  It  has  been  hypothesized  that  these  modes  should 
couple  to  and  thus  produce  resonant  absorption  of  irradiating 
electromagnetic  energy. 

(1)  Resonances:  Past  Results 

In  1984  and  1985  Edwards,  Davis,  Saffer  and  Swicord3 > * 
(EDSS)  reported  the  observation  of  relatively  sharp  peaks  in  the 
microwave  absorption  spectrum  of  aqueous  DNA  solutions  in  the 
frequency  region  of  roughly  1  to  10  GHz.  These  observations  were 
surprising  and  controversial  because  it  had  been  believed  that 
the  strong  coupling  between  the  polymer  and  the  viscous  solvent 
would  damp  any  resonant  modes.1*2  Attempts  to  replicate  these 
results  both  in  our  laboratory3  and  elsewhere6  have  led  to  the 
conclusion  that  they  are  experimental  artifacts,  unrelated  to  the 
normal  "organ-pipe"  modes  of  the  DNA  molecule.  It  seems  apparent 
that  the  above  mentioned  reservations  about  the  damping  effect  of 
solvent-DNA  coupling  on  the  existence  of  microwave  absorption 
peaks  were  well-founded. 


( 2 )  Microwave  Measurements 

The  existence  of  resonant  modes  in  polymers  has  been 
observed  in  crystalline  fatty  acids  and  indeed  in  normal  al¬ 
kanes.7  We  have  adopted  the  Dorf man-Van  Zandt  hypothesis  that 
the  failure  to  observe  these  modes  in  DNA  in  solution  stems  from 
their  being  quite  heavily  damped  through  strong  coupling  to  their 
viscous  surroundings.  The  important  quantity  in  determining  this 
coupling  is  the  product  of  the  mode  frequency  with  the  viscous 
relaxation  time  of  the  solvent.  For  aqueous  solvents,  this 
product  is  on  the  order  of  0.1,  and  thus  this  coupling  is  rather 
strong  and  as  a  result  the  vibrational  modes  are  not  observed. 

The  difficulty  is,  of  course  that  lower  viscosity  (and  thus 
shorter  relaxation  time)  solvents  are  simply  not  available. 
However  it  should  be  possible  to  reduce  the  coupling  between  the 
DNA  modes  and  the  solvent  by  increasing  the  viscosity  of  the 
latter.  In  the  oversimplified  but  still  representative  case  of  a 
single  viscous  relaxation  time,  the  damping  of  a  vibrational  mode 
is  proportional  to 
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1  +  (  COT)  “  . 

Observe  that  this  function  is  a  maximum  for  u>t  =  l .  Note 
particularly  that  for  wt  >>  1,  this  factor  falls  off  as  1/tuT; 
thus  by  increasing  the  solvent  viscosity,  n»  and  consequently  the 
relaxation  time  t  (since  n  =  G-t),  the  damping  can  be  decreased. 

We  have  initiated  a  program  to  investigate  the  effect  of  solvents 
whose  relaxation  times  can  be  varied  by  changing  the  temperature 
to  see  if  these  longitudinal  resonances  an  be  detected. 

( 3 )  Raman  Scattering  Studies 

The  same  kind  of  thinking  allows  one  to  recognize  that  if 
one  increases  u  sufficiently,  a  similar  result  can  be  achieved. 

To  accomplish  this,  we  have  also  initiated  a  program  to  study  the 
Raman  spectra  of  shorter  polymers  of  the  type  studied  by  Vergoten 
et  al .  such  as  stearic  acid.  In  contrast  to  their  work  which  was 
carried  out  in  the  pure  crystals,  our  studies  would  be  performed  in 
relatively  dilute  solutions.  The  aim  of  this  work  is  twofold:  (1) 
to  establish  that  the  damping  of  the  resonant  modes  does  indeed 
vary  with  frequency  and  relaxation  time  as  suggested  above,  and 
(2)  to  establish  the  appropriate  parameters  (e.g.,  t)  for  a  range 
of  solvents  as  a  function  of  such  properties  as  temperature  and 
salt  concentration. 
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A.  Microwave  irradiation  system  for  exposure  of  biological  cell 

at  2.45  GHz. 

1.  System  description 

The  original  exposure  system  described  in  section  1.2.1  of  the 
first  annual  report  was  modified  to  consist  solely  of  coaxial 
components  except  for  the  applicator.  The  term  'applicator'  refers 
to  the  section  of  waveguide  were  the  sample  is  loaded  for  irradia¬ 
tion.  A  block  diagram  of  the  system  is  shown  in  Fig.  1.  The 
incident  power  and  the  reflected  power  are  measured  at  the  outputs 
of  the  first  and  second  20  dB  coaxial  couplers  respectively.  The 
transmitted  power  is  measured  at  the  output  of  the  applicator. 
Since  the  sample  represents  a  mismatch  in  the  line,  optimum  energy 
coupling  is  obtained  by  adjusting  the  double  stub  tuner  to  minimize 
the  reflected  power. 


Fig.  1:  Block  diagram  of  system  for  microwave  exposure  of 

biological  cells  at  2.45  GHz. 


2. 


Applicator  design 

To  alleviate  some  problems  with  temperature  uniformity  and 
settling  of  the  cell  suspensions,  the  original  applicator  was 
modified  such  that  the  sample  tube  extends  out  through  both  the  top 
and  the  bottom  of  the  waveguide.  This  configuration  leads  to  a 
more  uniform  electric  field  distribution  within  the  portion  of  the 
sample  tube  inside  the  waveguide  and  allows  the  use  of  a  magnetic 
stirring  bar  which  rotates  at  the  bottom  of  the  sample  tube  to 
maintain  the  cells  in  suspension.  The  action  of  the  magnetic 
stirring  bar  also  reduces  temperatures  gradients  within  the  sample. 
Temperature  measurements  along  the  long  axis  of  the  sample  tube 
showed  temperature  uniformity  to  within  .2°C  at  a  nominal  SAR  of 
400  mW/g.  1  cm  O.D.,  3.6  ml  polyethylene  sample  tubes  with  screw 
cap  tops  are  used  as  sample  containers.  Cooling  is  provided  by 
circulating  silicone  oil  DC200/5  through  a  sealed  bath  surrounding 
the  sample  tube.  This  fluid  is  nontoxic,  chemically  inert,  and 
essentially  lossless.  Fig.  2  shows  a  longitudinal  cut  along  the 
center  of  the  applicator.  The  sample  temperature  is  monitored 
during  exposure  with  a  fiber  optic  Luxtron  probe  inserted  trough 
the  cap  of  the  sample  tube. 


Fig.  2:  Longitudinal  cut  along  applicator  loaded  with  sample 

tube . 
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Fig.  3:  Applicator  configuration  for  irradiation  of  attached 

cells . 

A  plexiglass  box  was  built  around  the  applicator  to  provide 
some  atmospheric  isolation  and  allow  flow  of  a  controlled  atmos¬ 
phere  of  95%  air  and  5%  COz  which  was  found  to  be  necessary  to 
achieve  pH  stability  during  long  time  runs,  e.g.  overnight  ex¬ 
posures.  For  greater  effectiveness,  the  air/C02  mixture  is 
delivered  into  a  polystyrene  vial  positioned  upside  down  over  the 
sample  tube  inserted  in  the  applicator. 

The  same  applicator  was  used  for  exposure  of  cell  monolayers 
attached  to  a  glass  substrate.  In  this  case  a  teflon  insert  is 
placed  in  the  sample  tube  to  hold  the  glass  slide  containing  the 
sample.  The  applicator  is  positioned  as  shown  in  Fig.  3.  The 
sample  temperature  is  monitored  during  exposure  with  a  fiber  optic 
Luxtron  probe  inserted  through  a  septum  in  the  cap  of  the  sample 
tube  and  located  directly  over  the  cell  monolayer. 

Applicators  identical  to  those  described  above  but  not 
connected  to  a  microwave  signal  were  used  for  sham  exposures. 
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SAR  measurements 


3  . 


The  Specific  Absorption  Rate  (SAR)  may  be  expressed  as 
SAR=CP  (dT/dt)  where  CP  is  the  specific  heat  of  the  irradiated 
sample  and  dT/dt  is  the  rate  of  increase  of  temperature  per  unit 
time  within  the  irradiated  sample.  Since  temperature  i-  .  .-asured 
as  a  bulk  property,  it  is  clear  that  the  SAR  determined  using  this 
definition  is  a  measure  of  energy  deposition  in  the  bulk  medium  and 
not  necessarily  in  the  cells  suspended  in  the  medium  or  the  cells 
bathed  by  the  medium  for  the  case  of  attached  cells.  Beaiing  this 
in  mind,  the  SAR  is  not  sufficient  to  specify  a  given  experimental 
condition  without  knowledge  of  the  volume  of  medium  exposed.  This 
is  particularly  true  for  the  case  of  exposure  of  monolayers. 

SAR  measurements  for  the  cell  suspension  case  were  carried  out 
by  monitoring  dT/dt  in  response  to  an  incident  signal  of  15  Watts 
at  various  positions  along  the  center  axis  of  the  sample  tube 
within  the  waveguide.  The  SAR  measured  at  points  where  the  sample 
tube  crosses  the  walls  of  the  waveguide  was  lower  by  about  40%  than 
that  measured  at  the  center.  This  is  probably  due  to  heat  loss 
since  the  ends  of  the  sample  tube  were  not  insulated  during  the 
measurement.  If  the  true  value  of  the  SAR  is  taken  to  be  the  value 
measured  close  to  the  center  of  the  sample  tube,  the  SAR  for  this 
configuration  can  be  approximately  expressed  as:  SAR  = 
Pn t t  [mW] / [g] .  This  represents  a  value  of  the  SAR  within  the 
waveguide.  The  average  SAR  would  be  lower  since  the  sample  extends 
outside  the  waveguide  and  circulates  within  it  at  a  rate  determined 
by  the  rotation  of  the  magnetic  stirrer.  The  average  SAR  is 
estimated  to  be  approximately  three  fourths  of  the  measured  SAR 
since  at  any  given  time  only  three  fourths  of  the  sample  are  inside 
the  waveguide . 

Similar  measurements  were  carried  out  for  the  cell  monolayer 
case.  The  sample  container  was  completely  filled  with  nutrient 
medium  making  sure  no  air  bubbles  were  left  inside.  Again,  the  SAR 
measured  at  points  where  the  sample  tube  crosses  the  walls  of  the 
waveguide  was  lower  by  about  40%  than  that  measured  at  the  center. 
The  SAR  for  this  configuration,  taken  to  be  the  value  measured  at 
the  center  of  the  tube,  can  be  approximately  expressed  as:  SAR  = 
1 . 2*Pn  e  t  [mW] / [g] . 

It  must  be  noted  that  SAR  values  for  experiments  with  cell 
suspensions  reported  in  the  quarterly  reports  covered  by  this 
reporting  period  were  specified  based  on  SAR  measurements  performed 
on  the  original  applicator  configuration.  With  the  new  applicator 
energy  coupling  into  the  sample  is  considerably  more  efficient. 
According  to  the  new  SAR  calibration  the  reported  SAR  values  should 
be  increase  by  a  factor  of  three.  Thus,  experiments  nominally 
carried  out  at  100  mW/g  were  actually  carried  out  at  300  mW/g. 
Similarly,  the  SAR  values  reported  for  experiments  with  attached 
cells  using  the  new  applicator  should  be  increased  by  a  factor  of 
4.8. 
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3.  Summary  of  studies  carried  out  using  the  2.45  GHz  irradiation 
system. 

Various  studies  were  carried  out  during  this  reporting  period 
using  the  2.45  GHz  irradiation  system  with  both  CW  and  pulsed 
microwaves.  These  are  summarized  below.  Details  of  each  experi¬ 
ment  are  reported  elsewhere  in  this  report. 

1.  Effects  of  2.45  GHz  CW  microwaves  on  two  interferon  regulated 
enzymes,  2' -5'  oligoadenylate  (2-5A)  synthetase  and  2-5A  dependent 
endoribonuclease  (RNase  L)  in  murine  L929  cells.  Exposures  were 
conducted  for  4  hours  at  37 °C. 

2.  Studies  of  CW  microwave  induced  changes  in  cell  cycle  profiles 
between  37°C  and  41°C. 

3.  Studies  of  microwave  induced  dome  formation  in  confluent 
monolayer  cultures  of  LLC-PK  pig  kidney  epithelirl  cells  with  both 
CW  and  pulsed  microwaves. 

4.  Effects  of  CW  and  pulsed  microwaves  on  cell  viability  and 
ornithin-decarboxylase  and  RNase  L  activities  in  HL-60  human 
leukemic  lymphocyte  cells  in  suspension.  Short  term  and  long  term 
(24  hr . ) 

exposures  were  conducted  at  37°C. 

C.  Studies  of  the  effects  of  low  level  microwave  radiation  on  the 
activity  of  enolase. 

A  series  of  experiments  were  carried  out  at  Howard  University 
in  collaboration  with  Dr.  S.  K.  Dutta  to  study  possible  effects  of 
low  level  (0.12  mW/g)  microwave  radiation  on  the  activity  of 
enolase,  an  important  enzyme  of  the  glycolithic  path.  This  work 
involved  the  use  of  a  Crawford  cell  with  compartments  within  which 
2  small  culture  flasks  (25  cm2)  could  be  placed,  one  at  either  side 
of  the  center  conductor.  The  experiments  were  carried  out  with 
NG108  clone  cells  grown  attached  to  one  of  the  large  surfaces  of 
the  flasks  and  used  when  a  confluent  culture  was  obtained.  The 
cells  bathed  in  5  ml  of  nutrient  medium  were  irradiated  in  the 
Crawford  cell  such  that  the  microwave  signal  passed  first  through 
the  cell  layer  and  then  through  the  medium.  All  experiments  were 
performed  at  915  MHz  with  16  Hz  (80%)  modulation  for  a  30  minute 
exposure  at  37°C.  Two  exposed  and  two  control  samples  were  used. 

After  irradiation  the  enolase  assay  was  conducted  on  sham  and 
exposed  samples  by  spinning  down  the  cells  in  a  table  centrifuge, 
homogenizing,  removing  the  cell  debris  by  centrifugation  and  using 
the  remaining  cell  extract  to  react  with  2-phospo-glyceric  acid  in 
a  reaction  buffer  to  form  phosphoenolpyruvate ,  a  compound  with  an 
absorption  peak  in  the  UV.  The  kinetics  of  the  reaction  were 
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observed  for  3  to  6  minutes  by  monitoring  with  a  spectrophotometer 
the  increase  in  the  absorption  peak  of  phosphoenolpyruvate  at  240 
nm.  To  obtain  the  specific  enolase  activity  in  (gM 
product/min) / ( gg  protein)  the  absorption  at  280  nm  corresponding  to 
the  protein  concentration  was  also  measured. 

Our  observations  indicated  that  several  corrections  and 
improvements  to  the  original  experimental  method  were  needed  in 
order  to  validate  the  results.  A  main  point  of  contention  was  the 
validity  of  the  cellular  protein  determination.  Since  the  cells 
were  maintained  in  a  medium  containing  approximately  10%  protein, 
we  advised  that  a  rinsing  step  was  desirable  prior  to  homogeniza¬ 
tion  to  remove  the  extracellular  proteins.  Rinsing  with  polyphos¬ 
phate  buffer  saline  by  mixing  with  a  vortex  was  found  to  be 
undesirable  due  to  possible  rupturing  of  the  cells  by  the  vortexing 
action.  Although  several  methods  were  tried,  it  was  not  possible 
to  devise  an  acceptable  technique  in  the  time  available  to  effec¬ 
tively  carry  out  this  step.  The  protein  assay  method  which 
consisted  of  measuring  the  absorption  at  280  nm,  also  came  under 
question  due  to  interferences  from  absorptions  present  in  that 
region  of  the  spectrum  due  to  other  entities  present  in  the  cell 
extract,  such  as  nucleotides.  It  was  suggested  that  the  Bradford 
method  for  protein  determination  be  used.  This  method  which 
utilizes  a  reaction  with  proteins  yielding  a  product  with  an 
absorption  at  595  nm  is  generally  considered  to  be  more  accurate. 

Thirteen  experiments  were  carried  out  five  of  them  without  the 
additional  rinsing  step.  For  the  latter  case  the  results  show  a 
consistent  increase  in  the  enolase  specific  activity  of  the  exposed 
samples  compared  to  the  controls  (Table  I) .  This  evidence  has  led 
us  to  conclude  that,  in  spite  of  experimental  inconsistencies,  the 
observed  effect  appears  to  be  real.  However,  it  should  be  pointed 
out  that  there  is  a  need  for  better  controlled  experiments  in  this 
area . 


TABLE  I. 


Ratio  of  the  specific  activity  of 
over  the  control  samples. 


enolase  in 


the  exposed 


Exp  ID 


E/C1 


E/C  2 


1 

2 

3 

4 

5 


2.245  ±  0.735 
1.233  ±  0.180 
1.440  ±  0.114 
1.448  ±  0.335 
1.122  ±  0.064 


1.304  ±  0.305 
1.102  ±  0.102 


1.  Protein  determined  by  absorption  9  280  nm. 

2.  Protein  determined  by  Bradford  method. 
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D.  Microwave  irradiation  system  for  exposure  of  cell  monolayers 

using  a  Crawford  cell. 

The  system  described  in  section  1.2  for  exposure  of  cell 
monolayers  can  be  used  to  irradiate  cell  cultures  of  up  to  2.5  cm2. 
Bioassay  techniques  used  to  investigate  the  effects  of  the  radia¬ 
tion  at  specific  end  points  often  require  a  larger  number  of  cells 
than  that  which  can  be  found  in  a  2.5  cm2  monolayer.  To  overcome 
this  difficulty,  a  new  microwave  irradiation  facility  was  set  up  to 
allow  exposure  of  cell  monolayers  of  up  to  100  cm2.  This  system 
makes  use  of  a  Crawford  cell  for  the  exposure  chamber.  The  system 
block  diagram  is  shown  in  Fig.  4.  With  the  components  presently  in 
use  the  system  can  operate  between  800  and  1000  MHz  with  a  maximum 
incident  power  of  10  Watts  limited  by  the  use  of  a  10-Watt  TWTA 
unit.  The  signal  generator  (HP  612A)  can  be  operated  in  CW  mode, 
with  internal  modulation  of  400  Hz  or  1000  Hz,  and  with  external 
modulation  including  pulsed  signals. 

The  Crawford  cell,  model  IFI  CC110,  was  fitted  with  shelves 
positioned  parallel  to  the  electric  field  at  either  side  of  the 
center  conductor.  Two  doors  for  sample  insertion  were  placed  on 
either  side  of  the  Crawford  cell  perpendicular  to  the  center 
conductor,  and  on  opposite  sides  of  it.  The  Crawford  cell  mounted 
vertically  on  a  rotating  stand  was  placed  in  an  incubator.  Rotary 
coaxial  joints  were  connected  at  both  input  and  output  ports  of  the 
Crawford  cell  to  allow  free  rotation  of  the  cell  permitting  access 
to  the  sample  insertion  doors  from  the  front  of  the  incubator. 
During  exposure,  sample  containers  are  placed  symmetrically  at 
either  side  of  the  center  conductor  to  maintain  symmetry  of  the 
electric  field  within  the  Crawford  cell. 


Fig.  4: 


Block  diagram  of  system  for  microwave  exposure  of  cell 
monolayers  using  a  Crawford  Cell. 
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SAR  measurements  were  carried  out  using  25  cm2  tissue  culture 
flasks.  Four  flasks  were  used  placing  two  at  either  side  of  the 
center  conductor  of  the  Crawford  cell.  The  effective  SAR  at 
various  power  levels  was  calculated  by  making  measurements  of  the 
incident,  reflected,  and  transmitted  power  to  compute  the  net  power 
loss  under  two  conditions:  a)  Crawford  cell  loaded  with  four  empty 
tissue  culture  flasks,  and  b)  Crawford  cell  loaded  with  four  tissue 
culture  flasks  filled  with  5  ml  of  nutrient  medium.  The  SAR  in  the 
medium  was  determined  from  the  ratio  of  the  difference  of  the  net 
power  loss  under  conditions  a)  and  b)  to  the  total  volume  of  the 
sample  (20  ml).  A  linear  regression  was  performed  on  the  measured 
data  to  obtain  the  following  functional  relationship  between  SAR 
(mW/g)  and  incident  power  (mM) :  SAR=3 . 74*10" 3 *Pi N c .  Results  are 
shown  in  Fig.  5.  It  must  be  noted  that  these  SAR  measurements  are 
valid  only  for  the  condition  described  above  and  are  representative 
of  power  deposition  within  the  bulk  nutrient  medium  and  not 
necessarily  within  the  cell  monolayer. 


SAR  vs. Pine, 4-25cm' flasks  f  5  ml  medium 


Fig.  5:  SAR  (mW/g)  vs.  Incident  power  (mW)  for  Crawford  Cell. 
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Temperature  control  was  provided  by  a  National  water  jacketed 
double  door  incubator.  Temperature  stability  within  this  incubator 
is  within  ±  0.1°C,  as  measured  with  a  Luxtron  fiber  optic  ther¬ 
mometer.  In  order  to  insure  that  a  stable  temperature  has  been 
reached,  samples  must  be  loaded  in  the  Crawford  cell  at  the  same 
temperature  as  the  incubator  temperature  (experimental  temperature) 
and  allowed  to  equilibrate  for  at  least  one  hour  prior  to  irradia¬ 
tion  . 


Additional  temperature  measurements  were  carried  out  to 
determine  the  operating  SAR  range  within  which  no  additional 
temperature  control  is  needed  to  maintain  temperature  stability. 
The  results  of  these  measurements  are  shown  in  Table  II.  These 
results  indicate  that  under  the  described  conditions  the  system  can 


be  operated 

at  SAR's 

below  about  5 

mW/g  without 

inducing  an 

appreciable 

temperature 

increase  within 

the  sample. 

TABLE 

II.  Temperature  increase  at 

various  SAR  levels. 

SAR 

Equil  temp 

Temp  increase 

in  10  mins. 

(mW/g) 

(°C) 

(°C) 

4.7 

38.65 

0.05 

8.3 

38.7 

0.30 

27.5 

41.3 

1.20 

E.  Test  fixtures  for  measurements  of  dielectric  response. 

Measurements  of  the  dielectric  response  of  essential  compounds 
present  within  a  biological  cell  such  as  DNA  and  other  polynucleo¬ 
tides,  offer  a  means  to  gaining  some  understanding  of  the  basis  of 
the  interaction  of  electromagnetic  fields  with  biological  matter. 

Two  test  fixtures  have  been  designed  to  perform  these  measure¬ 
ments,  one  for  frequencies  below  40  MHz  and  the  other  one  for 
frequencies  above  45  MHz.  A  Hewlett  Packard  4194A  Impedance/Gain 
Phase  Analyzer  was  used  for  measurements  in  the  range  between  100 
Hz  and  40  MHz.  A  major  problem  for  measurements  at  lower  frequen¬ 
cies  is  electrode  polarization,  particularly  when  measuring 
conducting  samples.  This  effect  can  sometimes  be  present  even  at 
the  low  MHz  range.  This  problem  can  be  corrected  for  by  using 
movable  electrodes.  Bearing  this  in  mind,  a  Hewlett  Packard  16451A 
dielectric  cell  with  movable  electrodes  and  specially  modified  to 
accept  liquid  samples  was  selected  to  make  the  low  frequency 
measurements.  Initial  measurements,  however,  have  shown  that 
additional  steps  are  necessary  to  overcome  the  electrode  polariza¬ 
tion  effect.  A  commonly  used  procedure  is  to  coat  the  electrode 
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surfaces  with  platinum  black.  While  this  cell  may  be  adequate  once 
the  electrode  surfaces  are  properly  treated,  it  requires  a  rela¬ 
tively  large  sample.  A  separate  cell  will  be  built  for  smaller 
samples  based  on  a  design  by  Dr.  Kenneth  Foster  (private  communica¬ 
tion)  . 


Beyond  a  few  hundred  MHz  dielectric  measurements  can  no  longer 
be  made  with  the  impedance  bridge  technique.  A  viable  approach  is 
the  use  of  coaxial  lines  which  can  generally  be  used  between  40  MHz 
and  18  GHz.  This  technique  makes  use  of  the  reflected  signal  from 
a  sample  located  at  a  defined  reference  plane  along  the  coaxial 
line.  A  wide  range  of  sample  configurations  are  possible.  In  all 
cases  the  sample  is  placed  at  the  end  of  a  coaxial  line  of  charac¬ 
teristic  impedance  Zo .  The  sample  may  constitute  part  of  the 
coaxial  line  or  be  placed  or  be  placed  at  the  end  of  the  line.  The 
latter  case  is  that  of  an  open  ended  coaxial  line.  For  our 
purposes,  a  dielectric  cell  was  designed  based  on  a  GR-900  connec¬ 
tor.  The  connector  design  was  modified  such  that  the  teflon 
support  for  the  inner  conductor  extended  past  the  end  of  this 
conductor  forming  a  cup  to  hold  the  sample.  The  outer  conductor 
which  was  flush  with  the  end  of  the  teflon  support  was  shorted  with 
a  cap  thus  closing  the  cup  to  seal  the  sample  and  complete  the 
circuit.  A  filling  chamber  connected  by  two  small  holes  with  the 
sample  chamber  was  built  in  the  shorting  cap  and  sealed  with  a 
rubber  septum.  This  allows  loading  of  the  sample  chamber  after 
creating  a  vacuum  to  eliminate  air  bubbles.  Measurements  with  this 
cell  were  conducted  with  a  Hewlett  Packard  8510  Network  Analyzer. 
Initial  results  indicated  that  the  practical  useful  range  is  below 
1.5  GHz.  An  open  ended  coaxial  line  will  be  used  to  probe  the 
dielectric  response  at  frequencies  above  1  GHz. 
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II.  BIOLOGICAL  APPROACH  TO  THE  PROBLEM 


Our  initial  research  efforts  on  effects  of  microwave 
irradiation  on  cultured  cells  involved  a  toxicological,  or  "cell 
injury"  approach.  This  approach  was  based  on  the  idea  that  events 
such  as  DNA  synthesis  in  S  phase,  mitosis,  and  survival  following 
plating  at  low  density  are  dependent  on  a  complex  assemblage  of 
metabolic  events,  and  that  perturbation  of  any  component  event 
would  thus  affect  the  observed  end  point.  These  efforts  were 
emphasized  in  Year  1,  and  are  summarized  in  that  annual  report. 
Additional  examinations  of  related  endpoints,  using  proliferating 
suspension-adapted  cells,  HL60  cells  and  a  transport  epithelium 
cell  line  are  summarized  in  this  report.  Perturbations  of  cell 
proliferation,  plating  efficiency  or  visibility  were  not  found  after 
microwave  exposure.  Examination  of  cell  cycle  timing  and  of 
membrane-related  differentiation  functions  also  showed  no  effects 
of  microwave  exposure.  These  observations  represent  a 
clarification  of  a  confusing  situation  in  the  literature,  where 
microwave  effects  have  been  reported,  but  without  adequate 
distinction  between  thermal  suid  athermal  conditions.  Given  the 
data  accumulated  to  date,  it  would  seem  that  athermal  effects  of 
microwave  irradiation  are  probably  subtle,  and  best  approached  at 
the  molecular  level. 

Two  specific  enzyme  systems  were  selected  for  analysis  during 
year  2.  The  rationale  for  selecting  these  systems  is  detailed 
below.  Additionally,  a  high  resolution,  two  dimensional  gel 
electrophoresis  system  was  developed  to  allow  examination  of 
synthesis  of  many  polypeptides.  The  selected  enzymes  and  the 
electrophoresis  systems  are  complementary,  in  that  both  specific 
endpoints  and  general  endpoints  are  being  approached  on  the 
molecular  level. 

1.  Pse  of  Proliferating  Suspension  Cultures,  and  Crawford  Cell 

Bffigflagfll. 

Previous  experiments  were  done  using  cells  which  had  been 
released  from  their  growth  surface,  placed  temporarily  into 
suspension  for  irradiation,  and  then  returned  to  monolayer  growth 
for  subsequent  assay.  These  cells,  which  had  not  been  adapted  to 
permanent  suspension  culture,  were  not  actively  proliferating 
during  the  temporary  suspension  conditions  employed  for  microwave 
exposure.  In  order  to  obtain  cells  that  would  actively  proliferate 
in  suspension  conditions,  L929  monolayer  cultures  were  placed  at 
high  density  into  spinner  flasks  with  low  calcium  medium,  and  cell 
density  was  allowed  to  decrease  as  adaptation  occurred,  and  active 
growth  was  observed. 
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Suspension-adapted  L929  cultures  were  exposed  to  conditions 
previously  used  for  irradiation  of  temporarily  suspended  cells. 
These  included  exposures  to  300  mW/g  average  SAR  microwaves  at 
37  *C,  as  well  as  exposures  at  temperatures  of  39,  41  and  42 *C.  In 
no  case  was  any  differential  effect  observed  in  cell  proliferation, 
plating  efficiency  or  viability  as  the  result  of  microwave 
exposures.  Actively  proliferating  cells  thus  were  not  more 
susceptible  to  microwave  effects  than  cells  temporarily  suspended. 

Installation  of  the  Crawford  Cell  by  Engineering  provided 
needed  facilities  for  irradiation  of  proliferating  monolayer 
cultures.  The  previously  used  wave  guide  devices  provided 
sufficient  space  for  use  of  suspension  cultures,  or  of  small 
coverslip  monolayer  cultures.  The  coverslip  cultures  provided  too 
few  cells  for  adequate  analysis  by  biochemical  or  flow  cytometry 
techniques.  Availability  of  the  Crawford  cell  allowed  exposures 
at  0.9  GHz  and  with  cell  quantities  sufficient  for  biochemistry  and 
cytometry . 

Initial  exposures  were  conducted  at  5  mW/g  SAR  to  confirm 
results  previously  obtained  with  suspension  cultures  at  2.45  GHz 
conditions.  Exposures  at  this  SAR  yielded  a  15%  increase  in  cell 
number  of  microwave  over  control  cultures  after  72  hrs.  As  a 
consequence  of  these  observations  the  engineering  group  made 
measurements,  and  determined  that  slight  elevation  of  temperature 
(less  than  0.5*  C)  was  occurring  at  5  SAR.  Subsequent  exposures  at 
SAR  of  1  mW/g  were  done,  conditions  under  which  no  heating  was 
measured,  and  showed  no  increase  in  cell  number  of  exposed  over 
control  samples.  Installation  of  a  fern  in  the  incubator  minimized 
temperature  buildup  within  the  Crawford  cell  due  to  microwave 
absorption,  but  any  subsequent  long  term  exposures  (greater  than 
approx .  4  hrs .  )  were  conducted  at  SARs  of  3  mW/g  or  less  to  avoid 
heating  artifacts. 

2.  Cell  Cycle  Analysis 

Methods 

Cultures  of  the  HL60  cell  line  were  inoculated  at  1x10s 
cells/mL  and  samples  were  taXen  for  counts  of  cell  density  and  cell 
cycle  analysis  at  intervals  of  twenty-four  hrs.  These  samples 
provided  baseline  information  needed  for  accurate  setup  of  log 
phase  cultures  for  irradiation  experiments.  Typical  DNA  profiles 
of  log  phase  cultures  are  represented  by  the  non- irradiated  samples 
as  shown  in  Fig.  1  and  Table  1.  Results  from  this  work  indicated 
that  cultures  established  24  hr.  prior  to  irradiation  would  double 
in  cell  density,  to  2xl05  cells/mL,  by  the  time  of  irradiation  and 
be  in  log  phase  growth  during  the  course  of  the  experiment. 

Cells  to  be  analyzed  by  flow  cytometry  were  processed  for 
propidium  iodide  ONA  staining  according  to  the  methods  outlined  in 
Protocol  #  27,  from  Ortho  Diagnostics.  This  protocol,  which 
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involves  lysis  of  the  cell  with  non-ionic  detergent,  was  selected 
since  it  allowed  immediate  processing  of  cells  without  the  need  for 
extended,  RNase  digestion,  and  since  it  minimized  cell  clumping. 
DNA  content  was  analyzed  for  10,000  cells  per  sample  using  an  Ortho 
Cytof luorograf .  Each  sample  was  analyzed  in  quadruplicate  for  sham 
and  irradiated  cultures,  and  in  duplicate  for  non- irradiated 
controls.  On  the  basis  of  the  DNA  profile,  regions  were  selected 
to  represent  Gl,  S  and  G2/M  phases  of  the  cell  cycle,  and  the 
percent  of  total  cells  in  each  stage  was  calculated. 

Irradiation  conditions  consisted  of  0.9  GHz  microwaves  which 
were  modulated  (80%)  at  60,  400  or  1,000  Hz,  with  an  SAR  of  2  mW/g. 
These  irradiations  were  done  for  24  hr.  in  a  Crawford  cell,  with 
cells  maintained  in  sealed,  25  cm2  flasks.  Sham  irradiated  cells 
were  placed  outside  the  Crawford  cell  on  a  shelf  in  the  incubator, 
and  non- irradiated  controls  were  maintained  in  the  C02  incubator  in 
which  the  cells  were  grown  prior  to  the  experiment. 

Results 

Typical  DNA  profiles  for  non-irradiated,  sham  and  microwave 
irradiated  samples  are  shown  in  Fig.  1.  Profiles  were  remarkably 
consistent  from  experiment  to  experiment,  and  within  experiments. 
Table  1  compares  the  cell  cycle  distributions  obtained  for  the 
three  irradiation  conditions  applied.  In  no  case  was  there  a 
significant  difference  within  an  experiment  to  indicate  an 
alteration  of  cell  cycle  progression  due  to  microwave  irradiation. 
Variation  for  a  given  cell  cycle  stage  between  sham  and  irradiated 
samples  was  within  1%  for  all  experiments,  and  these  samples 
displayed  no  stage  specific  variation  of  more  than  2.0%  compared 
to  the  non-irradiated  controls  which  were  maintained  in  a  separate 
incubator.  Cell  cycle  distribution  was  thus  unaffected  by 
microwave  exposure  over  a  24  hr.  interval.  Cell  viabilities  and 
densities,  determined  at  the  end  of  each  experiment,  were 
comparable  among  all  three  samples  for  each  experiment. 
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Table  1:  Cell  Cycle  Distributions  After  Microwave  Irradiation 


Modulat ion 

Saapla 

G1 

G2/H 

60  Hz 

HE 

11.6 

8.6 

79.8 

SH 

11.8  10.9 

8.4  10.379.4  *1.1 

HU 

1*2.2  10.4 

8.4  10.879.4 

ll.l 

400  Hz 

NE 

21.3 

12.6 

66.1 

SH 

19.6  10.4 

10.7  *0.2 

69.7  *0.4 

HU 

19.7  10.9 

10.9  10.5 

69.4  *1.0 

1000  Hz 

NE 

18.8 

7.4 

73.8 

SH 

18.5  10.7 

7.1  *0.274.4  *0.7 

HU 

17.7  10.8 

7.1  10.175.2  *0.8 

NE  ■  nonaxposad  control,  SH  *  thaa  axposad,  HU  ■  aicroaava  axposad 


3.  Cell  Differentiation  Experiments 

Porcine  cells  of  the  LLC-PK  line  represent  one  of  several 
established  mammalian  cultures  which  exhibit  differentiated 
characteristics  of  a  transport  epithelium  when  grown  to  confluency 
(11).  Confluent  monolayers  establish  apical  tight  junctions 
between  adjacent  cells,  and  transport  water  passively  across  the 
monolayer  as  the  result  of  active  transport  of  salts  and  sugars. 
Accumulation  of  water  beneath  the  monolayer  produces  a  build  up  of 
positive  hydrostatic  pressure,  which  manifests  itself  as  a 
loosening  of  local  areas  of  the  monolayer  from  the  growth  surface. 
This  produces  transient  "blisters"  or  "domes"  in  which  areas  of  the 
monolayer  elevate,  and  subsequently  lower  back  onto  the  growth 
surface,  presumably  as  hydrostatic  pressure  ebbs  with  leakage  of 
water  through  the  monolayer.  The  domes  can  be  counted  under  the 
microscope,  and  the  number  of  domes  per  microscope  field  used  to 
quantify  the  extent  of  transport  activity  at  a  given  time.  This 
situation  thus  allows  for  microscopic  assessment  of  a  complex  set 
of  coordinated,  differentiated  functions.  Agents,  such  as  the 
adenosine  analog  theophylline,  which  interfere  with  membrane 
transport  activity ,  thus  diminish  the  extent  of  dome  formation  ( 9 ) . 

The  transport  phenomena  exhibited  by  LLC-PK  cultures  are 
primarily  the  result  of  membrane  mediated  functions.  Since  at  least 
some  alterations  of  membrane  transport  have  been  associated  with 
microwave  exposure  ( 10 ) ,  it  was  felt  that  the  epithelial  culture 
system  would  provide  a  means  of  assaying  for  altered  membrane 
functions . 

Methods :  LLC-PK  cells  were  maintained  as  monolayer  cultures 
in  Eagle's  Minimum  Essential  Medium  supplemented  with  10%  fetal 
bovine  serum  and  20  mM  HEPES  buffer.  Stock  cultures  were  routinely 
split  prior  to  confluency,  so  as  to  avoi  problems  of  diminished 
cell  growth  rate  following  differentiation.  For  irradiation 
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experiments  cell  suspensions  were  seeded  at  lxlO6  cells  per  dish, 
into  petri  dishes  coverslips.  Coverslips  were  harvested  for 
experiments  after  1  or  2  days  growth,  depending  upon  whether 
confluent  or  subconfluent  specimens  were  desired.  Irradiations 
were  conducted  for  24  hr.  using  either  continuous  wave  conditions 
at  300  mW/g  average  SAR,  or  1  nsec  duration  pulses,  at  100  pulses 
per  second,  at  am  SAP.  of  25  mW/g.  Irradiations  were  conducted  in 
the  temperature  control  device  previously  described,  with  sham 
samples  held  in  an  identical  device  but  without  microwave  exposure. 
Control  cultures,  maintained  in  the  original  incubator,  were 
assayed  for  dome  formation  prior  to,  and  immediately  after,  each 
experiment.  Dome  formation  was  assessed  by  fixing  coverslip 
preparations  in  2%  glutar aldehyde  in  phosphate  buffered  saline  to 
preserve  the  dome  conditions  immediately  at  the  end  of  exposure, 
and  then  counting  domes  per  field,  using  phase  contrast  optics, 
for  an  area  of  40  mm2.  Variation  in  the  number  of  domes  per  area 
was  found  between  different  platings  of  cells,  but  was  consistent 
within  an  experiment.  To  assess  the  ability  of  the  cells  to  show 
response  to  perturbations  by  altering  dome  formation,  some  cultures 
were  treated  with  0.1  mM  theophylline. 

Results ;  Application  of  theophylline  to  confluent  LLC-PK 
cultures  reduced  dome  formation  by  40% ,  with  control  populations 
showing  an  average  of  45  domes  per  field,  and  theophylline  treated 
cells  27  domes  per  field. 

Irradiation  experiments  were  conducted  of  both  confluent  and 
subconfluent  cultures,  with  the  idea  that  irradiation  might 
interfere  with  either  establishment  of,  or  maintenance  of,  a 
transport  epithelium.  In  none  of  the  experiments  conducted,  with 
either  continuous  or  pulse  wave  exposures,  was  a  significant 
difference  in  dome  formation  observed.  Representative  values  are 
given  below  for  two  exposures. 

Table  2:  Dome  Analysis  of  LLCPK  Cultures 

Confluent  Culture.  Pulsed  Wave 

Average  Domes /Field  After  Exposure 


Sham 

18 

Exposed 

20 

Pre-Control 

29 

Post-Control 

22 
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Subconfluent  Culture.  Pulsed  Wave 


Average  Domes /Field  After  Exposure 

Sham  12 

Exposed  11 

Post-Control  22 

Cultures  maintained  in  the  original  incubator  consistently 
showed  more  domes/field  than  those  transferred  to  the  sham  or 
irradiation  chambers.  This  was  not  due  to  failure  of  cells  to  grow 
under  experimental  conditions  since  subconfluent  cultures  reached 
conf luency  in  the  irradiation  chambers ,  and  initiated  dome 
formation.  The  constant  flow  of  medium  across  the  coverslips  in 
the  chambers  amd  different  conditions  for  gas  exchange  may  explain 
the  difference,  as  may  the  transferring  of  cultures  required  for 
insertion  in  the  chambers.  Standard  deviation  for  control  cultures 
ranged  from  ±3-6  domes  per  field.  In  no  case  was  the  difference 
in  dome  density  between  sham  and  microwave  exposed  samples 
sufficient  to  indicate  an  effect. 

4.  Eng vme  Activities  and  Microwave  Exposures 

Ornithine  decarboxylase  (ODC),  am  enzyme  involved  in  the 
production  of  the  polyamines  putrescine  amd  spermidine,  has 
attracted  considerable  research  attention  (12).  ODC  is  one  of  the 
most  highly  inducible  amd  regulated  eucaryotic  enzymes,  and  is 
involved  in  pathways  importamt  to  RNA,  DNA  amd  protein  synthesis. 
ODC  activity  has  been  shown  to  vary  in  response  to  cell 
trams  formation,  and  correspondingly,  increases  with  cell  growth  amd 
proliferation.  Given  the  possibility  of  alterations  in  membrame 
electrical  activity  by  exposure  to  electromagnetic  fields,  amd  the 
fact  that  ODC  activity  responds  to  stimuli  via  membrame-mediated 
second  messengers,  it  has  been  selected  as  a  subject  for 
investigations  of  electromagnetic  field  effects.  Two  to  five  fold 
elevations  in  ODC  specific  activity  have  been  reported  for 
mammalian  cultures  exposed  to  60  HZ  electromagnetic  fields,  10 
mV/cm  for  1  hr  (1).  Increased  activity  has  also  been  reported  for 
exposure  of  cultures  to  450  MHz  microwave  irradiation  modulated  at 
10-20  Hz  (2). 

We  selected  ODC  as  am  enzyme  for  investigation  since  it  has 
been  demonstrated  to  respond  to  electromagnetic  effects,  and  due 
to  its  involvement  in  several  essential  biochemical  pathways. 

Additional  enzymes  were  selected  for  investigation  which  had 
not  been  previously  examined  for  responses  to  electromagnetic 
fields,  but  whose  importamce  to  the  functioning  of  the  immune 
system,  amd  whose  regulation  via  membrame  receptors,  suggested 
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fruitful  possibilities  for  investigation.  The  enzymes  selected, 
2 '  ,  5 ' -oligoadenylate  synthetase  (2-5A  synthetase),  and  the  2-5A- 
dependent  RNase  (RNase  L)  are  synthesized  by  cells  following 
response  to  interferon.  Upon  interferon  stimulation,  basal  level 
activities  are  elevated  from  three  to  several  thousand  fold  for  2- 
5-A  synthetase  and  up  to  30  fold  for  RNase  L,  depending  upon  the 
cell  line  and  other  factors  (15).  Elevation  in  2-5A  synthetase 
activity  results  in  synthesis  of  2'  to  5*  linked  oligoadenylates 
from  ATP.  This  2-5A  activates  RNase  L,  which  cleaves  cellular  or 
viral  single-stranded  RNA  at  sites  following  UU  or  UA  sequences. 
Approximately  24  hrs  post  stimulation  is  required  for  maximal 
activity  of  either  enzyme  to  be  realized  after  exposure  to 
interferon.  The  two  enzymes  have  been  implicated  in  the  anti-viral 
and  growth  inhibition  effects  of  interferon,  and  may  also  be 
involved  in  cellular  differentiation  and  growth  regulation  (7). 

ODC  Assays 

Methods :  Stock  cultures  of  L929  cells  were  grown  as 

monolayers  in  Eagle ' s  Minimum  Essential  Medium  supplemented  with 
10%  fetal  bovine  serum  and  20  mM  HEPES  buffer.  3.5  mL  of  cell 
suspension,  2xl06  cells/mL,  from  midlog  phase  cultures  were  used 
for  each  sham  or  microwave  irradiation.  Microwave  and  sham 
exposures  were  carried  out  at  300  mW/g  average  SAR,  37  *C,  for  4  hr 
in  the  wave  guide  apparatus  already  described.  Cells  were  then 
harvested  for  assay.  Cells  were  harvested  directly  from  monolayer 
for  the  "culture  phase"  ODC  assay. 

ODC  Specific  Activity  Measurements:  Cells  were  washed  3  times 
in  cold  PBS  and  then  frozen  at  -70*  C  for  storage.  Thawed  pellets 
were  suspended  in  lysis  buffer  (25mM  Tris,  pH  7.5,  2.5  mM 
dithiothreitol ,  0.1  mM  EDTA,  0.1%  Nonidet  P-40),  and  S-10  fractions 
prepared  by  centrifugation  at  10,000  rpm  for  15  min,  4’C.  Protein 
determinations  were  done  by  use  of  the  Bio-Rad  commassie  blue 
method.  100  /zg  of  S-10  sample  protein  was  used  in  a  250  uL  assay 
mixture  including  40  /xM  L-ornithine,  275,000  dpm  L-[UC] -ornithine 
(approx.  50  mCi/mMole),  4  uM  pyridoxal  phosphate,  1.25  mM 
dithiothreitol  and  5.0  mM  Tris-HCL,  pH  7.5.  The  reaction  was 
carried  out  by  incubation,  37 *C,  in  stoppered  test  tubes  carrying 
polypropylene  wells  containing  200  fiL,  1M  hyamine  hydroxide.  After 
60  min  the  reaction  was  stopped  by  injection  of  300  /xL  20%  TCA 
through  the  rubber  stopper.  Acid-released  HC02  was  dissolved  in 
the  hyamine  hydroxide  during  a  further  30  min  incubation  at  37 *C. 
The  well  containing  the  hyamine  hydroxide  was  removed,  the 
hydroxide  neutralized  with  5  fxL  glacial  acetic  acid,  and  the  well 
placed  into  a  scintillation  vile  for  counting  (13). 

Results :  ODC  activities  have  been  shown  to  vary  with  the 

state  of  cell  proliferation.  To  assure  consistent  conditions  for 
microwave  experiments  cultures  were  established  and  allowed  to 
progress,  over  five  days,  through  initial  lag  phase,  into  log  phase 
growth  and  into  stationary  phase.  ODC  activities  were  measured  at 
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0  and  4  hrs.  and  then  at  24  intervals.  Results  are  presented  in 
Table  3.  These  data  provided  baseline  ODC  values  for  L929  cultures, 
and  demonstrated  that  ODC  regulation  was  occurring.  Subsequent 
cultures  for  irradiation  were  accordingly  allowed  to  progress  to 
late  log  phase  for  use,  so  as  to  provide  cells  for  exposure  that 
exhibited  baseline  levels  of  ODC,  but  which  were  actively 
proliferating. 

Also,  to  be  sure  that  alterations  in  ODC  activities  by 
exogenous  agents  would  occur  with  the  L929  cell  line,  cultures  were 
exposed  to  phorbol  methyl  ester  in  concentrations  ranging  from  160 
to  1280  nM.  Regulation  was  observed,  with  the  640  nM  yielding 
approximately  3  times  higher  ODC  activity  than  comparable  controls. 


Table  3;  Variation  in  00C  Activities  with  Culture  Conditions 
Ti—  (hrs)  Phase  of  Growth  00C  Activity  (units/o  prof  in) 


0 

Inoculta  (Late  Log  Phasa) 

3.9  t 

4 

Lag  Phasa 

122 

24 

Early  Log  Phasa 

60.6  t 

48 

Mid  Log  Phasa 

13.3  t  5.7 

72 

Lata  Log  Phase 

3.7  i  1.0 

% 

Early  Stationary  Phasa 

3.1  i  1.1 

120 

Lata  Stationary  Phasa 

0.6  t  0.6 

14 

On*  unit  equals  the  Mount  of  enzyae  releasing  1  nanoaole  of  CO^  par  Minute 
at  37*C. 


Results  from  ODC  assays  of  3  preliminary  experiments  involving 
4  hr  sham  or  microwave  exposures  at  300  mW/g  SAR  were  initially 
compiled.  Specific  activities  of  microwaved  samples  were 
virtually  identical  to  those  of  sham  irradiated  samples,  with 
values  of  56  ±  17,  and  54  ±18  units/g  protein,  respectively.  As 
follow  up  to  these  data  a  series  of  exposures  of  L929  monolayers 
was  begun,  using  the  Crawford  cell  and  SARs  of  3-5  mW/g,  over  a 
time  course  of  0  to  8  hrs.  Additional  experiments  were  begun  using 
the  human  cell  line  HL60,  and  repeating  the  four  hour  exposures  in 
suspension,  using  continuous  wave  or  pulsed  microwaves.  These 
experimental  series  were  initiated  at  the  end  of  the  grant  year. 


b.  Interferon  ..System 

Methods:  Mid  log  phase  L929  monolayer  cultures,  grown  in  150 
cm2  flasks,  were  trypsinized  and  placed  into  suspension  at  a 
density  of  5xlOb  cells/ML.  Cells  for  microwave  or  sham  exposures 
were  placed  into  3.6  mL  tubes,  inserted  into  the  temperature 
controlled  waveguide  devices  previously  described,  and  maintained 
in  suspension  with  stirring  during  the  period  of  irradiation. 
Monolayer  cultures,  maintained  in  the  original  incubator,  were 
retained  for  subsequent  basal  level  control  values,  or  were  treated 
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with  Interferon  (100  IU  murine  alpha  plus  beta)  for  positive 
controls.  Irradiations  were  at  300  mW/g  average  SAR,  2.45  GHz,  for 
4  hours  at  37  *C.  Following  irradiation  cells  were  returned  to 
tissue  culture  flasks  and  incubated  at  37 *C  for  18  hr.  In  other 
experiments  cells  were  maintained  as  log  phase  monolayers  in  25  cm2 
flasks,  and  irradiated  at  5  mW/g  SAR  in  the  Crawford  cell  for  4  hr, 
and  incubated  an  additional  18  hr.  Cells  were  then  harvested, 
washed  in  cold  PBS  and  pellets  frozen  at  -70 *C.  Pellets  were  lysed 
and  S-10  fractions  prepared  as  in  the  ODC  assay. 

2-5A  synthetase  activity  was  assayed  by  isolating  2-5A 
synthetase  from  S-10  fractions  onto  poly(I)-poly(C)-cellulose, 
followed  by  incubation  with  4.0  mM  ATP  for  2  hr  at  37*C  (17).  2- 

5A  synthesis  was  measured  by  direct  competition  radiobinding  assay 
using  [32P]-2,5A  (6).  Units  of  synthesis  were  nMoles  of  2-5A 

synthesized/mg  protein/hr. 

RNase  L  activity  levels  were  monitored  by  binding  of  [^P]- 
2 , 5A .  100  /J.G  S-10  protein  was  incubated  with  5,000  cpm  [32 

P ] -2 , 5 A  overnight  at  0*C  (6).  Samples  were  filtered  through 

nitrocellulose  and  binding  measured  by  scintillation  counting  of 
filter-bound  radioactivity.  In  order  to  measure  RNase  L  specific 
activity  the  enzyme  was  immobilized  and  partially  purified  from  S- 
10  preparations  using  2-5A  bound  to  cellulose.  The  degradation  of 
poly (UJ-t^PjCp  (10,000  cpm/time  point  per  specimen)  was  assayed  in 
the  presence  and  absence  of  100  nM  trimer  2-5A  as  a  means  of 
demonstrating  enzyme  specific  activity.  20  mL  aliquots  were 
precipitated,  and  results  expressed  as  percent  of  original  cpm 
remaining  for  each  time  point  (14). 

Results:  Comparison  of  2-5A  synthetase  specific  activities 

from  microwaved  and  sham  populations  of  suspension  cells  and  from 
control  monolayer  cultures  revealed  no  significant  changes  in 
enzyme  activity  as  either  a  result  of  microwave  exposure,  or  of 
cell  handling  conditions  (Fig.  3).  Interferon  treated  monolayers, 
however,  exhibited  an  anticipated  increase  (approx.  191  fold)  in 
2-5A  synthetase  activity  from  163  to  31,200  units. 

Binding  of  2-5A  by  RNase  L  was  affected  both  by  conditions  of 
cell  handling,  and  by  microwave  exposure.  This  is  demonstrated  in 
Fig.  4.  Expressed  as  a  percent  of  control  binding,  the  sham- 
irradiated  cells  showed  an  approximately  19%  increase,  microwave 
irradiated  cells  an  approximately  42%  increase  and  interferon- 
treated  monolayers  approx,  a  67%  increase  in  binding.  Removal  of 
cells  from  the  growth  surface  and  maintenance  under  sham  conditions 
for  4  hrs  thus  produced  an  increase  in  RNase  L  binding  activity. 
Addition  of  microwave  exposure  more  them  doubled  this  binding 
activity,  but  did  not  reach  the  levels  achieved  when  interferon  was 
added  to  monolayer  cultures  subjected  to  neither  sham  nor  microwave 
conditions . 

Increased  RNase  L  activity  in  response  to  microwave  exposure 
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was  further  investigated  by  measuring  RNase  L  specific  activity  to 
assess  whether  increased  binding  actually  indicated  increased 
enzyme  catalytic  activity.  Results  of  these  measurements  are 
presented  in  Fig.  5.  Results  paralleled  those  for  2-5A  binding, 
with  highest  activity  displayed  by  interferon  treated  monolayer 
controls,  lowest  activity  by  non-treated  monolayers,  and  sham  and 
microwave  irradiated  samples  showing  respective  increases  over  non- 
treated  monolayers.  This  established  the  fact  that  the  more 
readily  performed  2-5A  binding  assay  for  RNase  L  could  be  used  as 
an  indicator  of  increased  enzyme  specific  activity. 

These  results  were  of  particular  interest  due  to  the  fact  that 
only  one  of  the  two  enzymes,  RNase  L,  displayed  increased  activity 
in  response  to  microwave  exposure.  Under  conditions  of  interferon 
stimulation  both  enzymes  would  show  increased  activity.  The 
increased  RNase  L  activity  of  the  sham- irradiated  cells  may  have 
resulted  from  their  removal  from  monolayer  and  4  hr.  of  suspension 
culture,  during  which  time  no  cell  proliferation  would  be  expected 
to  occur.  RNase  L  activities  may  be  altered  if  cells  are 
maintained  in  non-proliferating  conditions  (5,8).  Given  these 
facts  it  was  decided  to  irradiate  cells  maintained  in  monolayer 
culture.  Irradiations  for  4  hrs  at  5  mW/g  SAR  were  conducted  in 
the  Crawford  cell  and  RNase  L  binding  activities  were  measured 
after  18  hr  additional  incubation.  Results  from  six  experiments 
are  compiled  in  Table  4.  As  can  be  seen  from  the  table  the  results 
obtained  were  equivocal.  Average  values  were  slightly  higher  for 
the  microwaved  samples,  but  sample  to  sample  variation  was  high, 
giving  standard  deviations  too  large  for  any  conclusions  to  be 
drawn.  Lower  than  anticipated  protein  values  for  the  S-10  samples 
of  the  least  3  experiments  suggested  that  at  least  some  of  the 
variation  might  have  stemmed  from  unanticipated  problems  with  cell 
cultures  or  early  processing  of  the  cell  pellets.  At  the  end  of 
the  grant  year  these  experiments  were  being  repeated  to  minimize 
sample  to  sample  variation  and  so  determine  whether  the  effect 
could  be  duplicated  with  a  proliferating  population. 


Tatola  4:  RNata  L  Binding  Activitiat  of  Irradiatad  Honolavara 
Expri— nt  <  Parcant  of  Control  CPU  l^P)-2-5A  Bound  Haan  tSD 


Shaa 


1 

2 

3 

4 

5 

6 


128.0 

106.0 

71.1 

71.1 

54.1 
60.7 


81.8i29 


MW  1  137.0 

2  140.0 

3  139.0 

4  57.8 

5  54.1 

6  66.7 


99. It 44 
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5 .  2-D  Gel  Electrophoresis 


Given  the  difficulties  involved  in  hypothesizing  which  cell 
products  would  most  likely  be  affected  by  exposure  to  microwave  or 
other  electromagnetic  radiation,  a  more  global  approach  to 
examining  molecular  perturbations  of  the  cell  was  sought  to 
complement  the  work  on  specific  enzymes.  Two  dimensional  gel 
electrophoresis  ( 2DGE)  allows  high  resolution  of  polypeptides  and, 
combined  with  isotope  labeling  techniques,  provides  a  well  defined 
means  of  analyzing  many  translation  products  simultaneously  (Cellis 
and  Bravo,  1984). 

Facilities  and  techniques  for  performing  translation  assays 
by  2DGE  were  established.  Using  a  pH  3-10  gradient  for  the 
isoelectric  focusing  first  dimension,  and  10%  polyacrylamide  SDS 
slab  gels  for  the  second  dimension,  a  system  was  established 
capable  of  resolving  approximately  1100  different  HL60  polypeptides 
when  used  in  conjunction  with  silver  staining.  Labeling  of  cells 
with  35S-methionine  at  20-30  /xCi  per  3xl06  cells  in  methionine  free 
medium  provided  sufficient  specific  activity  in  30  to  60  min  to 
load  5x10s  cpm  or  more  per  gel.  This  activity  allowed  exposures  of 
2  to  4  days  for  adequately  resolving  approximately  600  polypeptides 
on  X  ray  film,  provided  gels  were  impregnated  with  a  water  soluble 
fluor  before  drying. 

Initial  experiments  involving  exposure  to  microwaves  were  done 
while  the  conditions  for  labeling  and  fluor ography  were  being 
developed.  HL60  cells  grown  in  25  cm2  flasks  were  used  as  controls 
or  exposed  samples  for  SAR  3  mW/g  continuous  wave  or  60  Hz  sine 
wave  modulated,  0.9  GHz  microwaves  for  24  hr.  Proteins  were 
stained  with  silver,  and  control  and  experimental  pairs  were 
examined  manually  to  see  if  spot  patterns  of  the  polypeptides  had 
changed.  No  significant  changes  were  found  in  these  samples,  but 
since  the  silver  technique  stains  all  proteins,  and  is  not 
selective  for  those  synthesized  during  exposure,  only  very 
pronounced  alterations  in  synthesis  might  be  detected.  For 
example,  the  new  appearance  of  a  large  quantity  of  one  or  a  few 
polypeptides  would  be  readily  observed.  Variations  in  synthesis 
over  shorter  time  intervals  might  also  not  be  detected.  Subsequent 
exposures,  initiated  at  the  beginning  of  year  3  have  employed 
isotopic  labeling  and  f luorography ,  so  that  only  the  polypeptides 
synthesized  during  exposure  will  appear. 

6.  Discussion 


Results  from  Year  2  verified  that  microwave  exposed  actively 
proliferating  cells  displayed  the  same  resistance  to  perturbations 
of  cell  growth,  viability  and  plating  efficiency,  and  the  same 
thermal  susceptibilities,  that  were  seen  with  monolayer  cultures 
temporarily  suspended  for  irradiation.  Further,  membrane  related 
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differentiation  functions  of  the  cell  line  LLC-PK  were  not  altered 
by  either  continuous  wave  or  pulses  microwave  exposures  over 
periods  of  24  hr. 

Crawford  Cell  irradiations  of  actively  proliferating  cultures 
of  the  HL60  cell  line  showed  no  alterations  of  cell  cycle 
distribution  between  sham  irradiated  specimens  and  those  exposed 
to  either  continuous  wave  or  modulated  microwaves  at  0.9  GHz. 
These  results  resolve  the  earlier  question  as  to  whether  cell  cycle 
shifts  were  observed  in  cultures  which  subsequently  displayed  very 
low  level  bacterial  contamination.  Healthy,  actively  growing  cells 
were  not  altered  in  cell  cycle  progression  by  24  hr.  of  microwave 
exposure . 

Of  the  two  enzyme  systems  selected  for  analysis  one,  the 
interferon  related  system,  showed  response  to  microwave  exposure. 
RNase  L  binding  and  2-5A  activities  showed  an  approximate  doubling 
over  sham  values,  and  approx.  40%  increase  overall  when  compared 
to  monolayer  controls .  The  sham  values  ( 20%  over  controls )  may 
result  from  cessation  of  cell  proliferation  during  the  4  hr 
irradiation  period.  The  question  that  arises  under  these 
circumstances  is  whether  microwave  irradiation  is  influencing  RNase 
L  gene  regulation  in  some  direct  way,  or  is  affecting  it  indirectly 
through  some  effect  related  to  cell  proliferation.  None  the  less, 
it  is  interesting  that  the  effect  observed  was  approx,  two  thirds 
that  seen  in  interferon  treated  monolayers ,  and  that  the  RNase  L 
increase  was  not  concomitant  with  a  corresponding  increase  in  2-5A 
synthetase  activity.  Additional  exposures  of  monolayer  cultures 
at  0.9  GHz  will  be  conducted  in  year  2  to  see  if  irradiation  of 
actively  proliferating  monolayer  cells  will  produce  the  same 
results,  a  possibility  suggested  by  the  initial  data  reported  here. 
Additional  assays  will  be  done  for  the  HL60  human  cell  line. 

The  regulation  of  ODC  activity  in  L929  cultures  was 
established,  and  the  information  used  to  initiate  assays  of 
microwave  exposed  samples.  Initial  results  showed  no  significant 
alterations  as  the  result  of  microwave  irradiation,  but  this  work 
will  be  continued  to  explore  a  number  of  exposure  conditions,  and 
to  investigate  results  in  both  the  L929  and  HL60  cell  lines.  The 
regulatory  flexibility  of  ODC  (12)  makes  it  am  attractive  candidate 
target  for  electromagnetic  exposures,  amd  will  be  further  pursued. 

The  two  dimensional  gel  electrophoresis  system  yielded  no 
detectable  differences  in  protein  patterns  when  silver  staining  was 
used  to  compare  sham  and  exposed  cultures .  Only  gross 
perturbations  would  likely  have  been  detected  under  these 
circumstances,  however,  and  isotopic  labeling  was  required  for  more 
sensitive  assay.  Labeling  amd  fluorographic  conditions  using  S-35 
methionine  have  been  established,  allowing  up  to  lxlO6  cpm  to  be 
loaded  onto  each  gel  for  detection  of  even  small  amounts  of 
proteins.  Experiments  were  begun  to  analyze  isotopically  labeled 
proteins ,  examining  results  initially  from  ELF  exposures  at  60  Hz , 
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and  then  progressing  into  modulated  microwave  exposures. 

With  the  desire  to  extend  analysis  further  into  the  molecular 
realm,  experiments  were  planned  for  Year  3  to  examine  general  RNA 
alteration  as  well  as  specific  mRNA  variation  due  to  ELF  and  to 
microwave  exposures.  The  experiments  planned,  as  well  as 
extensions  of  those  already  in  progress,  would  thus  cover  general 
and  specific  transcriptional  patterns,  general  translational 
patterns,  and  specific  enzyme  activities. 
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.gure  1:  DMA  distribution*  for  log  phase  HL60  cultures,  as 
determined  by  flow  cytometry.  The  T  axis  is  cell  number 
for  a  given  level  of  DMA  fluorescence,  and  the  X  axis  the 
relative  amount  of  ONA  as  determined  by  propidium  ioidide 
fluorescence.  Dotted  vertical  lines  indicate  the 
approximate  divisions  for  Cl.  s  and  G2/M  amounts  of  DHA 
Controls,  P«-irradiation;  lb.  Controls,  post- 

irIa?iatlCa;  lc*  Shaa  irr*diated;  Id.  Microwave  exposed 
0.9  GHz,  400  Hz,  80%  modulation.  24  hr. 


2-5A  SYNTHETASE  ACTIVITY 
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Figure  2.  2-5A  Synthesis  Activities  at  37°C.  Unit*  are  nanomoles 

of  2-5A  synthesized  per  ng  protein  per  hour. 
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2-5A  BINDING  TO  RNASE  L 
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Figure  3.  Binding  of  ,aP  labeled  2-5A.  Values  are  given  as 
percent  of  binding  by  control  aonolayer  cells. 
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Figure  4.  Degrade  cion  of  Poly  (tf)-{,2PlCp  by  BNase  L.  Symbol 
designations  are  as  follows:  circles  •  monolayer 
control;  upright  triangles  *  sham  irradiated;  squares  * 
microwave* irradiated;  and  inverted  triangles  »  interferon* 
treated  monolayers.  Closed  symbols  represent  assays  done 
in  the  presence  of  unlabeled  triaer  2-5A  (100  nM);  open 
symbols  represent  assays  done  without  2-5A. 
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